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ABSTRACT 


M. 


This  report  is  concerned  with  the  analysis  and  computer  simulation 
<3f  a tVi-moored  ‘oceanic  buoy  structure  having  "a  near -horizontal  cable 
connected  between  two  of  the  three  main^see  Figurel^legs . The 
Method  of  Imaginary  Reactions  is  used  in  conjunction  with  the 
Method  of  Successive  Approximationsto  determine  the  equilibrium 
configuration  of  the  cable  array.  A special  technique  is  developed 
to  extend  the  above  methods  to  include  cable  arrays  with  one  internally 
redundant  loop.  Computer  search  routines  are  developed  to  ensure 
convergence  to  the  solution  for  equilibrium  positions  of  the  cable 
array. 

The  analysis  permits  the  inclusion  of  discrete  elements  and 
floating  devices  distributed  along  all  cables.  Expressions  are 
derived  in  detail  for  the  hydrodynamic  forces.  Both  tangential  and 
normal  drags  are  included  in  the  calculation  of  these  forces;  however, 
the  forces  induced  by  wave  action  are  not  considered. 

A computer  program  for  implementing  the  analysis  is  included, 
and  a sample  output  appears  in  Appendix  I.  ^ 

The  work  reported  herein  was  conducted  during  the  period  of 
.September  1968  through  December  1970. 
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SYMBOLS  AND  NOTATIONS 


The  symbols  and  notations  used  in  this  report  are  defined  as  they 

appear  in  the  context.  The  most  important  ones  are  listed  here  for  . . . 

. . • « • •*  * * * 

, . • • *••••’  * 

‘reference : 


A.  Symbols  as  used  for  Tri-Moored 
( ®n » » cn ) 

m,n 

XTEN(m,n) 

Cc;.x  Cc;y  Cc;z 
m,n'  m,n'  m,n 

ce;x  ce;y 

k,m,n'  k,m,n 

CD 

^k,m,n 

c" 

m,n 

(? 

m,n 

»n 

n 

fm,n 

<„•  Fm,n> 


Structure 

the  coordinates  of  the  nth 
cable  anchor 

the  effective  cross-sectional 
area  of  the  (k,m,n)th  elemental 
device 

the  extensional  rigidity  of  the 
(m,n,)th  cable  segment 

the  drag  constants  of  the  (m,n)th 
cable  segment 

the  drag  constants  of  the 
(k,m,n)th  elemental  device 

the  coefficient  of  drag  of  the 
(k,m,n)th  elemental  device 

the  coefficient  of  drag  of  the 
(m,n)th  cable  segment  when  this 
segment  is  normal  to  the  stream 

the  coefficient  of  drag  of  the 
(M,N)th  cable  segment  when  this 
segment  is  parallel  to  the  stream 

the  diameter  of  the  (m,n)th  cable 
segment 

the  component  of  drag  force  per 
unit  length  in  the  (m,n)  direction 

the  components  of  the  external 
force  acting  at  the  (m,n)th 
cable  station 


the  components  of  the  imaginary 
reactions  applied  to  the  M(n)th 
stations  of  cables  2 and  3 


the  components  of  the  hydrody- 
namic force  per  unit  length 
acting  on  the  (m,n)th  cable  segment 


the  components  of  the  hydrody^ 
paudc.^fprce  acting  on  .the^ 
(k,m,n)th  elemental  device 


the  components  of  the  lumped 
drag  force  at  the  (m,n)th  cable 
station  due  to  the  distributed 
hydrodynamic  forces  along  the 
cables 


the  components  of  the  lumped 
drag  force  at  the  (m,n)th  cable 
station  due  to  the  hydrodynamic 
forces  on  the  elemental  devices 

unit  vectors  in  the  (x,y,z) 
directions,  respectively 


the  total  number  of  elemental 
devices  attached  to  the  (m,n)th 
cable  half-segment  adjoining 
the  (m-l,n)th  cable  station 


the  total  number  of  elemental 
devices  attached  to  the  (m,n)th 
cable  segment 


the  stressed  length  of  the 
(m,n)th  cable  segment 


the  index  of  the  mth  station 
or  segment  on  the  nth  cable 


the  total  number  of  stations  or 
segments  on  the  nth  cable 


(RX  , Ry  . Rz  ) 
m,n’  m,n  m,n 


the  ratio  of  drag  coefficients 
for  the  (m,n)th  cable  segment  = 

cp  /cN 

m,n  m,n 

the  components  of  the  resultant 
force  in  the  (m,n)th  cable  segment 


k,m,n 


the  stressed  distance  of  the 
(k',m,n)th  elemental  device  from 
the  (m-l,n)th  cable  station 


k,m,n 


the  unstressed  distance  of  the 
(k,m,n)th  elemental  device  from 
the  (m-l,n)th  station 


V(m,n  : £ , 


the  tension  in  the  (m,n)th 
cable  segment 

2 

the  integral  of  V (z)  along  the 
(m,n)th  cable  segment  from  £ 
to  E2  equal  to 

f\ V2  [Z(m,n:£)]d£ 

the  weight  (or  buoyancy)  per 
unit  length  in  water  of  the 
(m,n)th  cable  segment 


k,m,n 


the  weight  (or  buoyancy)  in 
water  of  the  (k,m,n)th  elemen- 
tal device 


U,y,z) 


the  lumped  weight  (or  buoyancy 
force  at  the  (m,n)th  cable  station 

fixed  Cartesian  coordinates 


(x  ,y  _ ) 

m,n  m,n  m^n 


the  coordinates  of  the  (m,n)th 
cable  station 


z(m,n:£) 


(a  *Y_  „) 

m,n  m,n  m,n 


the  parametric  representation 

of  z along  the  (m,n)th  cable 

segment  = _ 

Z . +v  _ 
m-l,n  'm,n£ 

the  direction  cosines  of  the 
(m,n)th  cable  segment 

the  sine  of  the  angle  between 
the  (m,n)th  cable  segment  and 
the  stream 


xii 


■ *m».0  .-air  i iirtwOwrilMilMii 


(AF*.  \ .AfL  \ .AF?,,  \ ) 

M(n),n’  M(n),n’  M(n),n 


WI,m,n 


B.  Symbols  As  Used  for  Tie  Lea 


XXTEN(m) 


cc:x,  cc:y,  cc:z 

m m m 

ce:x  ce:y,  ce:z 

m * m ’ m 


(»;,  f;,  kz) 

m m m 


( c:x  hc:y  hc:Z) 
m ’ m * m 


the  components  of  the  additive 
forces  applied  to  the  M(n)th 
stations  of  cables  2 and  3 

the  hydrodynamic  constant  of 
the  (m,n)th  cable  segment  = N 

PCm,n  m,n/2 

the  hydrodynamic  constant  of 
the  (k,m,n)th  elemental 
device  = pcD  A 

k,m,n  k,m,n/2 


the  effective  cross-sectional 
area  of  the  (k,m)th  elemental 
device 

the  extensional  rigidity  of  the 
m th  cable  segment 

the  drag  constants  of  the  m th  cable 

the  drag  constants  of  the 
(k,m)th  elemental  device 

the  coefficient  of  drag  of  the 
(k,m)th  elemental  device 

the  coefficient  of  drag  of  the 
m th  cable  segment  when  this 
segment  is  normal  to  the  stream 

the  coefficient  of  drag  of  the 
m th  cable  segment  when  this 
segment  is  parallel  to  the 
stream 

the  diameter  of  the  m th  cable 
segment 

the  component  of  drag  force  per 
unit  length  in  the  m direction 

the  components  of  the  external 
force  acting  at  the  m th  cable 
station 

the  components  of  the  hydrodynamic 
force  per  unit  length  acting  on  the  mth 
cable  segment 


,,e:x  ,e:y  ,e:z» 
(h  , h h ) 

m m m 


(HC:X,  HC:y,  HC:Z) 
m ’ m ’ m 


(S*!I.  H*‘y,  <:z) 


i.  J»  k 

(J,m) 

j(m) 

j(m) 

B(m) 

L(M) 

(m) 

D 

r 

m 


Sv 

k,m 


the  components  of  the  hydrodynamic 
force  acting  on  the  (k,m)th  elemental 
device 

the  components  of  the  lumped  drag  force 
at  the  m th  cable  station  due  to  the 
distributed  hydrodynamic  forces  along  the 
cables 

the  components  of  the  lumped  drag  force 
at  the  mth  cable  station  due  to  the 
hydrodynamic  forces  on  the  elemental 
devices 

unit  vectors  in  the  (x,y,z)  directions, 
respectively 

the  index  of  the  Jth  elemental  device 
on  the  mth  segment 

the  total  number  of  elemental  devices 
attached  to  the  mth  cable  half-segment 
adjoining  the  (m-l)th  cable  station 

the  total  number  of  elemental  devices 
attached  to  the  mth  cable  segment 

the  stressed  length  of  the  mth  cable 
segment 

the  unstressed  length  of  the  m th 
cable  segment 

the  index  of  the  mth  station 

the  ratio  of  drag  coefficients  for  tne 

mth  cable  segment  = c / c 

in  m 

the  components  of  the  resultant 
force  in  the  mth  cable  segment 

the  stressed  distance  of  the  (k,m)th 
elemental  device  from  the  (m-l)  cable 
station 

the  unstressed  distance  of  the  (k,m)tn 
elemental  device  from  the  (m-l)th 
station 
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ST 

m 


W 


c 

m 


W 


e 

k,m 


the  tension  in  the  mth  cable  segment 
2 

the  integral  of  V (z)  along  the  mth 
cable  segment  from  5^  to  E^ 

fZ*  V2[Z(m  : C (]  dLC 

the  weight  (or  buoyancy)  per  unit 
length  in  water  of  the  mth  cable  segment 

the  weight  (or  buoyancy)  in  water  of  the 
(k,m)th  elemental  device 


the  lumped  weight  (or  buoyancy:  force  at 
the  mth  cable  station 


U,y,z) 


fixed  Cartesian  coordinates 


the  coordinates  of  the  ® th  cable  station 


Z(m.O 


(“.b.yJ 

m ,n  m 


A 

m 


the  parametric  representation  of  z along 

the  mth  cable  segment  = Z . + v £ 

m-l  m 

the  direction  cosines  of  the  mth  cable 
segment 

the  sine  of  the  angle  between  the  mtn 
cable  segment  and  the  stream 

the  hydrodynamic  constant  of  the  mth 

cable  segment  = PC  d . 

m m/  2 

the  hydrodynamic  constant  of  the 
(k,m)th  elemental  device  = P D 

V A,  , /**N 

x,m  k,m/^ 

a parameter  defining  distance  along  the 
mth  cable  segment 


.»  n „) 


respectively,  unit  vectors  along  the 
mth  cable  segment , normal  to  both  the 
mth  cable  segment  and  the  stream,  and 
normal  to  the  mth  cable  segment  but  in 
the  plane  that  includes  this  segment  and 
the  stream 


C.  Symbols  that  axe  Common  to  feoth  the  Main  and  the  Tie  Leg  Arrays 


COMPD 


COMPE 


a cut  off  value  that  defines  the 
acceptable  completion  of  the  successive 
approximation  iteration 


a cut  off  value  that  defines  the 
acceptable  completion  of  the  Imaginary 
Reaction  iteration 


a positive  definite  error  function 


the  current  vector 


the  current  magnitude  at  a height  z 
above  the  bottom 


a positive  convergence  factor  navmg 
the  dimensions  of  force 


the  density  of  surrounding  fluid 


the  angular  dir  of  the  current  with 
respect  to  x axis 
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INTRODUCTION 


] 


The  problem  of  determining  undersea  cable  configurations  due  t applied 
loadings  at  known  positions  along  the  cable  has  received  considerable  attention 
in  the  past  (2, 5. 8, 9)*.  However,  the  solutions  have  pertained  generally 
to  particular  loading  conditions.  A general  closed-form  analytical  method  eolv-  < 
large  variety  of  complex  cable  systems  has  not  been  available,  primarily 
because  of  the  nonlinear  characteristic  of  the  differential  equations  describing 
these  systems. 

Alekseev ^ ^ dealt  with  the  problem  of  a free~ended  cable  from  a continuum 
point  of  view.  He  obtained  a three-dimensional  solution  to  the  equilibrium 
equations,  which  included  gravity  effects,  along  with  arbitrarily  applied 
forces  along  the  cable.  Pode^  also  used  the  continuum  approach  to  deed, 
with  the  above  problem.  However,  in  treating  a towed  body,  he  only  considered 
a special  case  of  the  general  problem.  Both  Alekseev  and  Pode  assumed  an 
inextensible  cable  so  that  the  exact  integration  of  the  equilibrium  equations 
could  be  obtained.  Therefore,  these  solutions  are  independent  of  the  materials, 
used  for  the  towing  cables. 

Walton  and  Polachek^*^  approached  the  problem  of  a free-ended  cable 
from  a lumped  parameter  point  of  view.  They  developed  a two-dimensional 
numerical  solution  for  an  inelastic  line  in  water  when  the  end  conditions 
are  assumed.  Paquette  and  Henderson^  followed  a procedure  similar  to  that 
of  Walton  and  Polachek  using  analog  computer  techniques.  The  cable  was 

^Numbers  within  parentheses  refer  to  references  given  on  page  i':-.. 
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considered  to  be  elastic  but  constraints  were  placed  on  the  equilibrium 

position  of  the  cable  stations.  Both  steel  and  nylon  cables  were  considered 

(2) 

O'Brien  considered  the  case  of  fixed  end  elastic  cables,  such  as 

transmission  lines.  An  exact  continuum  solution  to  the  problem  was 

obtained;  however,  the  forces  applied  to  the  cable  were  assumed  to  be 

constant  over  sections  of  prescribed  length  so  that  the  shape  of  the  cable 

segments  could  be  expressed  as  an  elastic  catenary. 

/ o \ 

Skop  and  O'Hara'  developed  a technique  called  the  Method  of 
Imaginary  Jteactions,  which  is  an  extension  of  classical  consistent 
deformation^ theory  to  a nonlinear  problem.  As  reported  in  Keference- 8 
& 9 the  method  applies  to  elastic  non-redundant  cable  systems  and  uses 
lumped  parameter  representations  of  the  external  forces.  The  technique 
employs  a numerical  analysis  method  that  uses  a set  of  straight  segments 
to  represent  the  cable.  This  enables  the  non-linear  system  to  be 
represented  by  a set  of  linear  equations.  By  prescribing  a simple  method 
of  varying  the  redundant  reactions,  an  iteration  technique  is  used  which 
converges  to  the  correct  reactions  (and  consequently  the  correct  static 
equilibrium  configuration).  The  method  is  globally  convergent  and 
converges  to  actual  reactions  from  any  set  of  initially  estimated  reactions. 
The  only  restriction  reported  by  the  authors  (Reference  8 & 9)  is  that  the 
method  is  not  applicable  if  internal  loops  (a  redundant  structure)'  exist 
in  the  cable  system. 

(7) 

Savage  and  Sniffin  have  analyzed  a tri-moored  subsurface  float 
with  neutrally  buoyant  mooring,  legs  (fee  Fig.  lX,  utilizing  a three- 
dimensional  solution  that  assumes  that  the  neutrally  buoyant  legs  have  a 
catenary  shape  under  the  effect  of  current.  They  also  assume  that  if  the 
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Method  of  Imaginary  Reactions  to  handle  the  redundant  structural  cable 


array.'  ' This  report,  summarizing  the  research  work  of  Pahuja  published 
in  1970,  describes  the  use  of  the  Method  of  Imaginary  Reactions  in  the 
analysis  of  a single  loop  cable  array.  The  appendix  contains  the  computer 
program  developed  which  has  been  compiled  on  an  IBM  360A0  facility. 


CHAPTER  I 


METHOD  OF  IMAGINARY  REACTIONS 

A.  Background 

A classical  method  for  the  analysis  of  indeterminate  structures 

is  the  method  of  consistent  deformation,  the  method  of  consistent  distortion  or 

. (U) 

displacement  as  developed  by  James  Clark  Maxwell , This  method  of 

indeterminate  reaction  analysis  utilizes  equations  of  compatibility 

of  the  structure  to  supplement  the  equations  of  equilibrium  to  obtain 

a solution  to  the  unknown  redundants.  The  following  assumptions  are 

made  in  this  type  of  method: 

1.  The  structure  is  assumed  to  be  linear,  i.e.,  the  loads  applied 
are  proportional  to  the  displacements  and  the  displacements 
are  relatively  small. 

2.  That  there  are  no  gross  structural  distortions  or  instabilities 
upon  the  release  of  redundant  reactions,  i.e.,  upon  the  release 
of  the  "redundant  reaction"  the  basic  geometry  of  the  structure 
should  not  change. 


The  above  two  conditions  are  generally  not  met  in  the  analysis  of  cable 
arrays. 

(8) 

The  Method  of  Imaginary  Reactions  overcomes  these  restrictions 
for  the  analysis  of  cable  systems.  It  is,  however,  a natural  extension 
of  classical  consistent  deformation  theory. 

The  method  uses  the  following  assumptions: 

1.  The  bending  stiffness  of  the  cable  is  neglected. 

4 

6 


2.  The  external  forces  acting  on  the  cable  arrays  are  "lumped1 


If  the  cable  length  is  divided  into  segments,  then  each  cable  segment 


between  stations  can  be  assumed  straight.  The  analysis  technique 


provides  the  equilibrium  configuration  of  the  system  (including  the 


effect  of  cable  stretching)  which  is  determined  uniquely  from  formulas 


that  are  functions  of  only  the  applied  forces  and  the  reactions.  This 


offers  a a advantage,  in  that  the  solution  is  now  dependent  upon  the 


number  of  external  redundant  reactions  and  not  on  the  number  of  stations 


which  are  arbitrary.  The  third  assumption  is  necessary,  because  even 


though  convergence  is  still  obtained,  when  this  condition  is  not  met 


the  configuration  is  no  longer  uniquely  described  by  the  method 


B.  Description  of  the  Method  of  Imaginary  Reactions 


This  method 


briefly  in  this  section.  Consider  a single  cable  array  of  two-dimensional 


applied  forces  (the  method  is  not  restricted  to  two  dimensions) as  shown 


■ 


In  Figure  3,  Rsx  and  Hsy  acting  at  point  S,  Frx  and  Fry  acting  at 
point  R,  Fqx  and  Fqy  acting  at  point  Q,  are  externally  applied  and  are 


known.  If  P be  designated  as  the  anchor  point,  then  the  reaction  components 
Rpx  and  Rpy  can  be  determined  simply  by  summation  of  forces  in  the  x and  y 
directions. 

Since  these  reactions  are  known,  the  tension  in  segment  PQ,  the  length 
of  segment  PQ  under  the  tension  and  the  coordinates  of  point  Q become  known 
in  that  order  by  the  use  of  elementary  statics.  Using  the  same  type  of 
analysis  the  coordinates  of  point  R and  S can  also  be  found.  Once  the 
position  of  these  points  is  known,  we  have  obtained  the  equilibrium 
configuration  of  the  free-ended  cable.  If  point  S happens  to  be  the  second 
anchor  point,  then  the  fixed-ended  cable  system  is  in  equilibrium  in  the 
desired  configuration.  In  general,  however,  we  are  not  at  the  desired  fixed 
end  anchor  point , pt . T , tut  rather  at  point  S . Therefore  anchor  point 
reactions  Rsx  and  Rsy  at  point  S are  not  the  actual  anchor  point  reactions  which 
place  the  cable  end  at  the  desired  anchor  point  (T).  The  Method  of  Imaginary 
Reactions  provides  the  methodology  for  finding  the  actual  reactions  at  point 
T which  will  force  point  S to  move  over  to  the  anchor  point  (T),  producing 
the  true  equilibrium  configuration  of  the  system. 

This  is  done  by  first  applying  a small  force  AR  at  point  S vector ially 
directed  in  a direction  from  S towards  T.  If  the  coordinates  of  point  T 
are  (Tx,  Ty)  and  those  of  point  S are  (Sx,  Sy),  then  the  ratio  of  the 
component  of  this  additional  force  [(AR)  x 8c  (AR)y]  is  given  by: 

(AR)x  _ Tx  - Sx 

TARTy  = Ty  - Sy  (4,1)  • 

*The  numbering  of  all  important  equations  is  based  upon  the  format 
(a,b)  where  a is  the  number  of  an  equation  in  a given  chapter  b. 

(Equations  (1,1)  - (3,1)  have  been  deleted  from  the  text.) 
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The  magnitude  of  these  components  is  defined  by: 

(AR)x  = (Tx  - Sx)  (5,1a) 

(AR)  = (Ty  - Sy)  (5,1b) 

2 2 

vhere  E = (Tx  - Sx)  + (Ty  - Sy)  is  the  measure  of  error  between  the 
equilibrium  position  calculated  using  an  estimated  or  guessed  reaction  and 
the  correct  equilibrium  position.  E is  always  positive  and  vanishes 
uniquely  only  when  the  correct  equilibrium  has  been  obtained.  In  other 
words,  E*0  only  when  Tx  = Sx  and  Ty  = Sy  simultaneously.  However,  to 
arrive  at  the  correct  equilibrium  configuration,  it  is  necessary  to  assure 
that  E will  become  vanishingly  small.  To  assure  convergence,  a positive 
number,  6 , having  the  dimensions  of  the  force,  is  introduced.  6 is  a 
convergence  parameter  used  to  select  the  magnitude  of  the  additive  force 
AR  (as  demonstrated  in  equations  (5,1))  in  a manner  that  the  cable  array 
will  converge  to  the  correct  equilibrium  configuration.  Referring  back  to 
equations  (5,1),  the  ratios  (Tx  - Sx)/Je  and  (Ty  - Sy)/*^E  are  of  bounded 
variation  (lying  between  -1  and  +l) , and  E will  become  vanishingly  small 
only  when  <5-K).  When  this  happens,  the  additive  forces  (AR)x  and  (AR)y  are 
zero,  and  we  have  the  correct  reactions  and  thus  the  correct  equilibrium 
configuration. 


The  algorithm  for  obtaining  a solution  is: 

1.  Make  a reasonable  engineering  guess  as  to  the  components  of  the 
reaction  at  the  redundant  anchor  (References  8,9,  & 11  discuss 
methods  for  making  these  guesses). 

2.  Release  the  redundant  anchor,  i.e., point  S,  while  maintaining  the  guessed 
component  reaction  forces.  This  results  in  the  creation  of  a 
free-ended  cable. 
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3.  Calculate  the  equilibrium  position  using  static  force  equilibrium 
and  compute  the  quantities  (Tx  - Sx),  (Ty  - Sy'  and  E. 

£ 

- U.  Choose  an  initial  value  for  6 so  that  the  value  of  (AR)x=  (Tx-Sx), 

and  (AR)y  = (Ty-Sy)  can  be  found.  6 can  be  chosen  initially 
& to  be  very  large,  since  it  will  of  necessity  become  smaller  as  the 

solution  proceeds  step  by  step. 

5.  Step  4 results  in  an  additive  force  equal  to  R'  = R + AR  acting  at 
the  assumed  free  end.  Apply  this  additive  force  and  find  the  new 
equilibrium  position  of  the  cable. 

6.  The  next  step  involves  comparing  the  value  of  E'  (the  new  measure  ox 
error)  with  E (the  previous  measure  of  error.).  Two  possible  situations 
can  exist: 

a.  If  E'<E,  then  a successful  step  has  been  made  because  the  aim  is 
to  reduce  the  new  value  of  E to  zero.  In  this  case  the  same  value 
of  6 is  retained  and  a new  value  of  the  additive  force  (AR)  is 
calculated  and  added  to  the  previous  R to  calculate  again  a 

new  equilibrium  position.  Repeat  the  process  until  a stage  where 
E'  > E. 

b.  If  E<E' , then  6 is  too  large.  The  new  values  that  resulted  in  the 
new  measure  of  error  E'  are  rejected.  Go  back  to  the  last  iteration 
which  gave  the  previous  measure  of  error  E.  Then,  reduce  3 by 
halving  its  value.  Then  proceed  in  a repetitive  manner  with 
steps  U through  6,  until  the  value  of  E is  arbitrarily  small, 

at  which  point  the  released  end  is  arbitrarily  close  to  the 

f 

real  anchor  position,  and  hence,  a solution  to  the  problem 
has  been  obtained.  Figure  U is  a flow  diagram  describing 
the  procedure. 
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Figure  4:  A Flow  Diagram  Describing  the  Method  of 
Imaginary  Reactions 
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(8) 

This  application  of  Imaginary  Reactions  is  globally  convergent  ; 

1. e.,  no  matter  what  the  initial  guessed  reactions  are,  the  calculation 
converges  to  the  correct  coordinates*. 

2.  LUMPED  PARAMETER  REPRESENTATION 

One  of  the  assumptions  made  during  the  development  of  the  Method 
of  Imaginary  Reactions  is  that  the  distributed  forces  that  act  on  the 
cable  array  can  be  represented  by  lumped  forces  at  specified  points 
on  the  cable.  This  assumption  is  required  if  the  paths  to  the  free  ends 
of  the  system  are  to  be  determined  by  elementary  statics  and  formulas  that 
express  the  elongation  of  the  cables  under  tension.  The  only  condition  that 
these  elongation  formulas  have  to  meet  is  that  the  cable  length  increases 
with  the  increase  in  tension  and  be  a single  valued  function;  therefore,  the 
elongation  formulas  can  be  either  linear  or  non-linear. 

Some  of  the  guidelines  for  a successful  lumped  parameter  representation 
core  presented  below: 

1.  Each  cable  in  the  array  is  represented  by  at  least  two  stations, 
one  at  each  end. 

2.  Each  point  of  discontinuity  in  a physical  property  of  a cable 
is  represented  by  a station.  Consequently,  each  cable  segment 
in  the  array  is  homogeneous  with  constant  characteristics . The 
cable  can  be  described  by  as  many  stations  (segments)  as  are 
necessary  to  represent  the  cable.  Segments  into  which  the  cable 
is  divided  do  not  have  to  be  of  equal  length  along  the  cable 
array. 

3.  The  method  of  analysis  has  no  direct  bearing  on  the  number  of 
stations  on  the  cable  array.  As  many  additional  stations  as 


*The  analytical  proof  can  be  found  in  Reference  8. 
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are  necessary  to  obtain  a satisfactory  approximation  to  the 
continuous  equilibrium  shape  of  the  array  are  used.  The  error 
function  E,  the  positive  6,  and  the  various  values  of  A R do  not 
depend  upon  the  number  of  stations. 

h.  The  method  by  which  the  external  forces  acting  on  the  array  are 
to  be  lumped  at  the  cable  station  is  completely  arbitrary.  In 
this  report  a half  segment  technique  suggested  by  Paquette  and 
Henderson ^ is  used. 

C.  APPLICATION  OF  THE  METHOD  OF  IMAGINARY  REACTIONS  TO  DETERMINE  EQUILIBRIUM 

CONFIGURATION  OF  A MULTI-MOORED  CABLE  SYSTEM 

In  this  section,  the  Method  of  Imaginary  Reactions  is  discussed  in 
summary  fashion  to  determine  the  static  equilibrium  configuration  of  a 
multi-moored  cable  system.  (See  Figure  1)  A detailed  discussion  of  the 
analysis  of  this  structure  is  contained  in  Reference  9 & 11  . A summary 

discussion  of  the  application  of  the  Skop  & O'Hara  method  follows  as  background 
to  the  extension  to  the  single  internal  loop  analysis. 

1.  Definitions 

A few  definitions,  which  are  helpful  in  understanding  the  further  text, 
are  listed  below. 

1.  Simple  cable  array 

A simple  cable  array  is  defined  as  a cable  system  which  contains 
no  closed  loops,  i.e.,  it  is  not  internally  redundant. 

2.  A branch  point 

A branch  point  of  a simple  cable  array  is  defined  as  a 
point  at  which  a single  cable  "splits"  into  two  or  more  cables. 

3.  Primary  and  Secondary  Anchors 

Each  cable,  in  any  cable  system,  is  denoted  by  a number 


(n).  The  cable  is  said  to  terminate  at  a branch  point.  The 


...  — — 


index  n=l  is  reserved  for  the  cable  that  is  attached  to  an  anchor 
point  in  the  system.  In  case  more  than  one  anchor  point  exists 

in  a system,  then  the  point  to  which  cable  number  1 is  attached 
is  called  the  Primary  Anchor  and  the  other  anchors  are  called 
Secondary  Anchors.  It  is  the  Secondary  Anchors  that  are 
released  and  subjected  to  imaginary  reactions. 

2,  Coordinate  System  and  Notation 

For  the  tri-moored  System  analysis,  a right-handed  (x,y,z)  cartesian 
coordinate  system  is  used.  The  z axis  is  defined  to  be  parallel  to  the 
direction  of  gravity. 

The  location  of  the  nth  cable  anchor  is  given  by  a(n),  b(n),  c(n). 
Each  cable  in  the  system  is  represented  by  M(n)  stations.  The  location 
of  the  mth  station  on  the  nth  cable  is  denoted  by 


^Xm,n,  Ym,n,  Zm,n,)  (21,1) 

where  m * 1,2  M(n). 

The  stations  are  counted  from  the  primary  anchor  to  the  branch  point 
of  the  array  along  cable  1,  and  from  branch  point  of  the  array  to  the 
secondary  anchors  along  cables  2 and  3 respectively  as  shown  in  Figure  5. 


Thus,  for  a branch  point 


X(1,2)=X(1,3)=X(M(1),1) 

Y(1,2)=Y(2,3)S:Y(M(1),1) 


Z(1,2)=Z(2,3)=Z(M(1),1)  (22,1c) 

and  coordinates  of  the  three  anchor  points,  both  primary  and  secondary 
are  denoted  by 

(X(1,1),  Y(l,l),  Z(l,l))  3 (a*b*c)  (23,1a) 


(22,1a) 

(22,1b) 


(23,1a) 

( ^ * ^2  * C2 ^ 

(23,1b) 

( a3 • b3,  c3) 

(23,1c) 

(X(M(3),3),  Y(M(3) ,3) , Z(M(3),3))  “ (a3’  V c3)  (23,1c) 

The  external  forces  acting  at  the  (m,n)  th  station  are  defined  by  the 

components  along  the  x,y  and  z axis: 

where  m = 1,2  M(n) 

and  n = 1,2,3 

Thus,  forces  acting  at  an  anchor  point  are  represented  by 

(F*M(n),n)’  *Ol(n),n)»  F(M(n),n))  (25>l) 

and  forces  at  a branch  point,  which  is  indexed  by  (M(l),l)  conventionally, 

are  denoted  by 

F(M(1),1)’  ^Md),!)’  F(M(1),1)  (26,1) 


3.  Static  Equilibrium  Configuration 

It  was  stated  earlier  that  for  a two-dimensional  cable  system  if  external 
forces  in  segments  are  known,  then  the  tensions  in  and  orientation  of  the 
cable  segments  can  be  uniquely  determined  by  ordinary  static  method.  This 
can  be  generalized  to  three  dimensions  and  multi-moored  cables.  The 
essence  of  this  statement  can  be  summarized  as  follows: 

1.  If  the  external  forces  applied  to  the  arrays  are  constant, 

X Y Z 

then  the  components  of  the  resultant  force  (R  m,n,  R m,n,  R m,n,) 
acting  in  the  (m,n)  th  cable  segments  can  be  given  in  terms 
of  the  applied  external  force  through  the  following  expressions, 
for  m = M(n)  and  n = 2,3 


x _ x 

M(n),n  M(n),n 

t>y  py 

M(n) ,n  = M(n) ,n 

and 

z z 

R M(n) ,n  = F M(n) ,n 

for  m = M(l)  and  n = 1,  the  branch  point. 


(27, la) 
(27,1b) 

(27,1c) 


& 

,Z 


£ 

,z 


^(1),!  = FM(1),1 


+ R 

+ R 


,x  + RX 

2,2  + R2,3 

(28,1a) 

Y Y 

2 ,2  + R2,3 

(28,1b) 

Z „Z 

2 ,2  + R2,3 

(28,1c) 

17 


I 


and  for  m = 2,3  M(n)-1  and  n = 1,2,3 


R = F + R 
m ,n  m ,n  m+1  ,n 


(29,1a) 


R*  - F*  + B?\  _ 
m,n  m,n  m+l,n 


(29,1b) 


RZ  = FZ  + RZ 
m,n  m,n  m+l,n 


(29,1c) 


Thus  starting  from  the  secondary  anchor  points  the  forces  are  very 

readily  summed  up  to  give  the  reactions  at  every  station. 

Once  the  resultant  force  has  been  found,  the  tension  T is  given  by: 

m,n 


Tm,n=/  (R*  )2  + (fJ  >2  + (R*  / 

m,n  m,n  m,n 


(30,1) 


The  tension  in  each  segment  results  in  the  elongation  of  segments 


and  the  stressed  length  is  obtained  by 


BL  =BL  f , . m ,n 

m»n  ° m»n  l1+  XTENlm ,n ) 


(31,1) 


where  BL  is  unstressed  length, 

o m,n 

The  next  step  is  to  find  coordinates  of  each  station.  To  start  the 
solution,  one  starts  from  the  primary  anchor  because  its  coordinates 
are  known  and  are  fixed.  The  coordinates,  then,  are  found  from  the 
following  equations: 


18 


of  cable  2 and  3 will  not  be  at  the  true  anchor  points.  As  a measure 


of  the  distance  of  these  ends  from  the  correct  anchor  points,  the 


positive  error  function  E is  defined  as 


E will  uniquely  vanish  when  the  true  anchor  point  has  been  reached 


” •— 

I 

where  primes  denote  the  new  imaginary  reactions  and  the  additive 
forces  are  defined  by 

AFM(n),n  " ^ an-XM(n),n] 

= i [VW>,n] 

z <5 

AFM(n),n  ^ ^cn  ZM(n),n^ 


(35,1a) 

i 

(35,1-b) 

(35,1c) 


The  quantity  6 as  defined  earlier,  is  the  positive  convergence  factor. 
Thus  the  correct  equilibrium  configuration  is  obtained  once  E -*■  0. 

Thus,  in  the  above  section  the  use  of  the  Method  of  Imaginary  Reactions 
to  find  equilibrium  configuration  of  a tri-moor ed  structure  has  been 
demonstrated.  From  here  the  method  can  be  extended  for  n number  of  cables 
easily. 

D.  APPLICATION  OF  METHOD  OF  IMAGINARY  REACTION  FOR  POSITION-DEPENDENT 
EXTERNAL  FORCES 


In  the  preceding  section  the  ability  of  the  Method  of  Imaginary 
Reactions  to  determine  the  equilibrium  configuration  of  a cable  array,  both 
single  and  multi-moored  has  been  demonstrated.  However,  the  use  of  this 
method  to  determine  equilibrium  configuration  by  simple  statics  was 
dependent  totally  upon  the  lumped  external  forces  being  constant. 

However,  underwater  cable  arrays  are  subjected  to  external  forces  due  to 
weight  and  buoyancy,  and  hydrodynamic  forces,  position-dependent,  like  drag 
forces.  They  depend  upon  both  the  orientation  and  depths  of  the  cable 
segments  and  on  depths  of  elemental  devices,  with  the  result  that  the  Method 
of  Imaginary  Reactions  is  not  entirely  applicable  to  the  cable  structure 
we  wish  to  deal  with  in  this  report. 
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CHAPTER  II 


A NEAR-HORIZONTAL  CABLE  ELEMENT  IN  A MULTI -MOORED  ARRAY 

The  major  thrust  of  this  investigation  has  been  to  study  and  analyze 

the  tri-moored  structure  that  has  two  of  the  cables  connected  by  a 

horizontal  cable.  (This  cable  is  called  a tie  leg  in  this  report). 

A technique  is  developed  in  this  report  that  enables  the  Method  of  Imaginary 

Reactions  to  be  used  for  such  a structure.  More  recently.  Shop  and 

O'Hara  have  extended  the  work  reported  herein  to  the  general  redundant 

(12) 

structure . 

A.  THE  BASIC  APPROACH 

Let  one  first  consider  the  two  legs  of  a tri-moored  buoy  system, 
to  which  a tie  leg  has  been  attached.  This  structure,  then,  looks  like 
an  ’'A" . The  "A"  structure  can  further  be  divided  into  two  coupled  parts: 
l)  the  lambda  (A)  structure  and,  2)  the  tie  leg  structure  . Neglecting 
the  bending  stiffness  of  the  cables,  there  are  only  two  conditions  of 
coupling  between  the  lambda  structure  and  tie  leg. 

1.  The  geometrical  compatibility  condition 

2.  The  force  balance  condition 

All  that  is  required  is  a simple  lambda  structure  subjected  to  the 
above  constraints  between  (A)  structure  and  tie  leg,  and  then  the  structure 
will  behave  like  an  "A"  structure. 

B.  THE  GEOMETRICAL  COMPATIBILITY  CONDITION 

Let  one  start  with  a lambda  (A)  structure  having  A as  an  apex  and 
B and  C as  anchor  points.  (See  Figure  6)  Points  D and  E represent 
approximately  the  position  at  which  one  wishes  a tie  leg  inserted  in  the 


are  the  total  resultant  forces  that  the 


structure . (FFY)at  D and  (FFY*)at  E 

lambda  structure  exerts  on  the  tie  leg.  Then  using  the  Method  of  Imaginary 

Reactions  the  equilibrium  configuration  of  the  structure  is  obtained  which, 

in  result,  gives  the  position  of  points  D and  E under  loads. 

If  points  D and  E are  known,  then  in  order  to  meet  the  geometrical 

compatibility  constraint,  the  tie  leg  chordal  length  must  fit  these  points. 

It  is,  therefore,  required  that  the  tie  leg  shall  have  these  points  as  end 

constraints.  The  Method  of  Imaginary  Reactions  is  then  applied  to  the  tie 

leg  by  releasing  the  E end  and  the  equilibrium  configuration  of  the  tie 

leg  is  obtained.  This  gives  one  certain  end  reactions  (EYl)^  D and 

(RY1*)  . „ which  the  tie  leg  exerts  on  the  basic  lambda  structure  at  points 
at  E 

D and  E . 


-(RYlL*  „ = (FFY). 


(RI1*>at  E * E 


C.  THE  FORCE  BALAHCE  COWDITIOB 

In  order  to  meet  the  second  constraint  of  force  balance,  it  is 
important  that  the  following  conditions  be  met: 

'at  D ‘ 'at  D (l,2a) 

(1,2b) 

When  this  condition  has  been  met , the  complete  solution  to  the  problem  has 
been  obtained,  i.e.,  the  lambda  structure  is  behaving  as  an  "A"  structure . 

However,  in  case  the  above  mentioned  condition  is  not  met,  then  the 
whole  process  must  be  repeated.  That  is  (FFY)  at  D and  (FFY*)  at  E is 
replaced  by  a new  assumed  reaction  (FFY)'  at  D and  (FFY*)?  at  E and  new 
position  of  points D and  E are  found  under  the  action  of  these  forces.  Then, 
these  positions  are  taken  as  reference  points  for  finding  the  equilibrium 
position  of  the  tie  leg,  which  gives  one  the  end  reaction.  This  process 
is  repeated  until  the  force  balance  condition  has  been  met.  Once  this 
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happens  the  two  constraints  have  been  taken  into  consideration  and  the 
problem  is  solved. 


" I j 

To  reduce  the  number  of  iterations  so  as  to  satisfy  the  above  constraint 
rather  quickly,  a convergence  algorithm  is  developed  in  Chapter  IV.  This 
algorithm  utilizes  the  values  of  FFY  and  RY1,  and  if  force  balance  constraint 

is  not  satisfied,  the  routine  employs  a technique  by  which  a new  value  of 

> > > > 

FFY'  (FFY  < FFY'  < RYl)  and  FFY*'  (FFY*  < FFY*'  < RY1*)  are  obtained.  This 

and  subsequent  values  of  FFY'  and  FFY*'  assure  convergence  to  the  correct  i 

values  of  forces,  which  satisfy  the  force  balance  condition,  in  minimum 

number  of  iterations. 


CHAPTER  III 


MODELING  AND  ANALYSIS  OF  THE  INTERNALLY  REDUNDANT  ARRAY  SYSTEM 

In  this  chapter,  a brief  description  of  the  techniques  used  to  model 
and  analyze  the  array  system  is  presented.  The  chapter  has  been  divided 
into  various  sections,  each  of  which  deals  with  one  phase  of  the  modeling 
process.  Further,  each  section  deals  separately  with  the  main  structure 
and  the  tie  leg  structure. 

A.  A FEW  GUIDELINES 

In  this  section  a few  guidelines  for  modeling  the  array  are  presented: 

1.  The  tri-moored  structure  consists  of  a subsurface  buoy,  anchored 
to  the  bottom  by  three  cables . These  cables  are  referred  to  by 
number  n where  n = 1,  2 and  3.  These  cables  are  broken  up  into 
M(n)  stations. 

2.  The  primary  anchor,  the  anchor  for  cable  n = 1,  is  represented 
by  (1,1).  All  other  stations  are  referred  to  by  subscript  pair 
(m,n)  where  m = 2,  3,  — M(n)  and  n = 1,  2,  or  3 except  for  the 
secondary  anchors  which  are  referred  to  by  (M(2),2)  and  (ml3),3). 

The  branch  point  at  the  subsurface  buoy,  in  this  case, is  represented 
by  (M(l),l).  This  is  the  same  point  as  (1,2)  or  (1,3)  but  by 
convention  is  represented  by  (M(l),l). 

3.  The  tie  leg  is  attached  to  cables  2 and  3,  at  stations  (L,2)  and  (L, 
respectively  where  L is  any  station  between  1 and  M(n).  The  tie 
leg  is  further  broken  up  into  MN  stations.  Both  arrays  are 
represented  diagrammatically  as  follows: 
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(1,0 


Figure  7 : Nomenclature  for  the  Tri-moored  Array  Structure  with  a Tie  Leg 

U.  Since  the  external  forces  acting  on  the  cables  have  to  be  lumped  at 
these  stations , each  cable  segment  between  the  stations  behaves  like 
a straight  line.  The  following  guidelines  are  used  for  lumping  these  forces 
i.  Each  point  of  discontinuity  in  a physical  property  of  a cable 
is  represented  by  a different  station.  Thus,  each  cable  in  the 
array  has  constant  physical  properties, 
ii.  As  many  additional  stations  are  used  as  necessary  to  obtain 

a successful  approximation  to  the  continuous  equilibrium  shape 
of  the  array.  However,  the  method  of  analysis  does  not 
depend  upon  the  number  stations  and  any  ultrafine  representa- 
tion of  the  cable  does  not  help. 

B.  DISCRETE  ELEMENT  NOMENCLATURE 

The  method  of  a ialysis  has  the  ability  to  take  into  account  discrete 
elements  in  the  cable.  These  discrete  elements  exert  a force  on  the  system. 

The  external  forces,  because  of  these  elements,  are  lumped  also  at  the 
stations  neighboring  them.  Half  segment  lumping  technique^  is  used  here. 

The  objects  which  are  attached  to  the  (m,n)th  cable  segment  on  the  main 

structure  are  indexed  by  (k,m,n)  where  k = 1,2  k(m,n),  counting 

ir.  the  direction  as  shown  in  Figure  8.  For  the  tie  leg — the  objects  aio 
indexed  by  (j ,m)  where  J = 1,  2 j(m). 
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By  convention,  the  device  on  the  station  itself  becomes  the  last 

object  on  the  segment.  Thus,  the  subsurface  buoy  becomes  the  last  device 

on  the  first  cable  and  is  indexed  by  (k(M(l),l)),  M(l),l). 

The  unstressed  distance  of  the  (k,m,n)th  device  from  the  (m,n)th 

station  is  given  by  S. 

k,m,n 

Similarly  in  the  case  of  a tie  leg,  the  element  (j*m)*th  would  be 

S.  ft,  away  from  the  mth  station. 

J ,m  ’ 

It  is  necessary  to  distinguish  between  forces  that  are  lumped  at 
the  lower  segment  from  the  forces  that  are  lumped  at  the  upper  segment. 
Half  segment  technique ^ is  used  to  lump  these  forces  on  the  respective 
stations. 

The  objects  that  are  attached  to  the  half  segment  adjoining  the 
(m,n)th  station  would  then  be  differentiated  from  the  objects  adjoining 
(m+l,n)th  segment  as  follows: 

Let  k(m,n)  represent  the  value  of  k such  that 

S,  <BL  /2 

k,m,n-  om,n 

for  k = 1,2  k(m,n) 

S.  >BL  /2 

k,m,n  om,n 

for  k = k(m,n)+l,  k(m,n) 


Similarly  for  a tie  leg  let  J (m)  represent  the  value  of  j such  that 


S.  <L/2 
J ,m  - m 


and 

for  J * 1,2 J(m) 
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(m-H.n ) 


(m  + l ) 


(k(m,n)+1,m,n) 


( j(m)+1,m) 


Cj(m)  ,m) 


(1,®) 


ELEMENTAL 

DEVICES 


(myn)  El«oents  on  the  main  arrays)  Elements  on  the  tie  leg 
Figure  8:  Representation  of  the  Discrete  Elemental  Devices 


These  equations  are  dlagrsamatically  represented  by  the  figures  shown  above 
C.  FORCES  ACTING  ON  THE  ARRAY  SYSTEM 


The  total  forces  acting  on  the  array  system  are  made  up  of  two  factors: 

1.  The  weight  and  buoyancy  forces — here  weight  and  buoyancy  of  discrete 
elements  as  well  as  of  the  cable  are  taken  into  account . 

2.  The  hydrodynamic  drag  forces  which  are  acting  because  of  the 
interface  between  the  cable  and  the  current  acting  id  water. 

1.  Weight  and  Buoyancy  Forces 

In  this  section  only  the  effect  of  weight  and  buoyancy  is  considered. 

The  weight  and  buoyancy  actually  can  be  differentiated  only  by  the  direction 
in  which  they  act.  The  weights  act  downwards  and  buoyancy  acts  upwards.  The 
Z direction  is  considered  positive  upwards. 

To  compute  the  weight  and  buoyancy  forces  let  the  weight  (or  buoyancy) 


per  unit  length  in  water  of  the  mth  cable  segment  be  given  by 


and  let  the  weight  buoyancy  in  water  of  the  (j,m)th  elemental  device  be 
given  as 


Then  to  lump  these  forces  as  a weight  (or  buoyancy)  force  W acting 

m 

at  the  mth  station,  the  half  segment  technique  is  employed.  That  is,  the 
distributed  and  discrete  forces  acting  on  the  half  segments  adjoining  the 


Technique 


Then  because  of  the  discrete  element  the  weights  that  are  lumped  at  station 
m,  acting  in  the  half  segment  of  the  (m-l)  segment  as  shown  in  Figure  9 are 
expressed  by 

J (m-l) 

j4(»-i)+l  3,m_1  (lj3) 

Similarly  the  forces  that  are  lumped  at  station  m because  of  weights 
acting  at  the  lower  half  segment  of  the  mth  segment  are  expressed  by 
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<— 1 1 OJ 


(2,3) 


j=j(m) 
£ W 


J=1 


e 

J »m 


Also,  the  lumped  force  at  station  mth  because  of  the  weight  ( buoyancy ) 
of  the  cable  would  similarly  be 


W . L . 
m-1  m-1 


+ 1 
2 “ 


m 


(3,3) 


Such  that  total  force  lumped  at  mth  station  due  to  weights  is  given  by 


„ _j(m-l) 
*“  £ 


m 


*?  . . + 1<?)  He 
J ,m-l  £ W, 

j=j  (m-l)+l  J=1 


+ 1 [ft  _L  . 
•g  m-1  m-1 


WCL  ] 
m m 


(4,3) 


Exactly  the  same  discussion  holds  good  for  lumping  the  weights  as  forces 
in  case  of  the  tri-moored  buoy  structure.  In  this  case,  however,  another 
factor  n,  the  number  of  cables,  has  also  to  be  taken  into  account. 

Then  for  the  tri-moored  structure 


w _k(m-l,n)  we  + k(m,n)  g 

m,n  , — ,£  . . . k,m-l,n  . £n  W 

k=k(m-l,n)+l  k=l  k,m,n 


+1  [WC  , BL  . + W BL 

t-  m-l,n  om-l,n  m,n  o m,n 


Where 

m = 2,3,4 M(n)-1. 

n = 1,2,3. 


(5,3) 
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and  for  m = M(l)  n = 1,  the  subsurface  buoy 

k(H(D,  1)  k(2,2) 


WM(1),1  = E wfc,M(l),l+  E wf,2,2 

k=k(M(l),l)+l  k=l 

k(2,3) 


k = 1 

♦ ![W£(1),1BLoM(1>-1  ^,2BLol.^i3BLoltj]  (6>3) 

which  are  various  expressions  for  representing  the  weights  as  lumped  forces 
on  the  stations.  In  deriving  this  expression,  however,  the  density  of  the  water 
has  been  assumed  to  be  constant. 

2.  Modeling  of  a Current  Profile 

Weight  and  buoyancy  forces  are  one  type  of  external  forces.  There  are 
additional  external  forces  which  act  on  the  cable  because  of  the  hydrodynamic 
interaction  between  the  cable  and  an  ocean  current.  This  interaction  produces 
drag  forces  which  are  dependent  upon  the  Velocity  and  direction  of  the  local 
ocean  current  i.e.,  the  angle  at  which  local  ocean  current  attacks  the  cable. 

The  method  of  analysis  developed  herein  makes  no  restrictions  on  the 
shape  of  the  current  profile.  For  the  sake  of  convenience,  however,  the  following 
assumptions  are  made: 

1.  The  ocean  current,  though  depth  dependent  in  magnitude,  is  uni-direct  i it  =»i 
and  normal  to  the  direction  of  gravity. 

2.  The  drag  force  component  which  acts  in  the  direction  normal  to  both 
the  stream  and  the  cable  is  zero. 

The  first  assumption  is  made  because  most  design  currents  are  given  as 
depth  dependent  and  uni-directional.  The  second  assumption  is  made  because 
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of  the  limited  experimental  data  available  to  support  a general  analytic 

• • • • ••  • • 

expression  for  the  side  component  of  "drag"  on  cables. 

The  first  assumption  dictates  that  the  current  possess  no  Z component.  If 
the  angular  direction  with  respect  to  the  x axis  is  denoted  by  4> , the 
expression  for  current  in  general  can  be  written  as 


V = V(Z)  [ i Cos  <fr  + J Sin  <J>] 


(7,3) 


Where 


V(Z)  = Magnitude  of  current  profile  at  a height  Z above  the  deepest  anchor  (Z=o), 


= Unit  vector  in  x direction 


J = Unit  vector  in  y direction 


3.  Hydrodynamic  Coordinates 


In  order  to  calculate  the  various  drag  forces  acting  on  the  system,  it  is 
convenient  to  transform  into  a new  set  of  natural  hydrodynamic  coordinates, 
defined  with  respect  to  direction  of  the  cable  and  the  current, 

The -following  analysis  is  presented  for  the  tri-moored  structure,  Howe- ec , 
the  same  approach  is  used  for  the  tie  leg  array. 

1.  In  order  to  find  a natural  hydrodynamic  coordinate  system,  we  proceed 
in  the  following  manner : 

i.  Let  n be  the  unit  tangent  to  the  (m,n)th  cable  segment, 
considered  positive  in  the  dir11  of  increasing  m 


ii .  Then 


V x xm,n 


(8,3) 


where  irm,n  is  a unit  vector  normal  to  both  the  (m,n)th  cable 
segment  and  the  current  and 

iii.  n m,n  = xm,n  x 7rm,n  is  a unit  vector  normal  to  the  cable  and  lying 
in  the  plane  that  includes  the  (m,n)th  cable  segment  and  the  stream 


In  order  to  express  the  hydrodynamic  base  reference  in  terms  of 
basic  x,  y,  z reference  frame,  we  obtain  the  following  relation: 


The  unit  vector  Tm,n  is  expressed  in  terms  of 


=a  i + 0 J+Y  k 

x m,n  m,n  m,n  u 'm,n 


(9,3) 


where  a m,n  0 m,n  and  Ym»Q  are  “the  cosines  of  the  (m,n)  th 


cable  segment  defined  by 
X -X 

nH-l.n  nun 


am,n=- 


BL 


m,n 

Y -Y 
m+l,n  m,n 


0m, n=- 


BL 


m,n 


(10,3a) 


(10,3b) 


Z , -Z 
m+l,n  m,n 


ym,n= 


BL 


m,n  (10, 3c) 

3.  To  express  mn,n  in  terms  of  basic  x,y,z  coordinates  we  proceed  as 
follows : 


nm,n= 


-*■ 

V x t m,n 
irm,n 


and 


Vxtm,n  = V(z)[i^  Cos  $ + J Sin  <fr]  x 

"t  t £ 
[om,n  ■*0m,nJ+Ym,n  ] 


=v(z) 


t 1 t 

Cos  $ Sin  $ 0 
om,n  0m  ,n  Y®,n 


■V(Z)[i  (ym,n  Sin  $)  -j”  (y^.p  Cos  4> ) 
+ k (0m,n  Cos  4>  -am,n  Sin  <|> ) ] 


3U 


Also 


$ X T 


m,n 


V 


T I A 
111,111  m,n 


Where 


Am  ^ is  the  sine  of  the  angle  between  the  (m,n)th  cable  segment  and  the 
stream  and  is  given  by 
A 


m,n"V  2 “7 

y + ( 


Ym  „ + (B  Cos  Sin  <J>)‘ 

m,n  m,n  m,n 


Also 


I = V(Z) 


And 


Thus 


m,n 


= /T 

n 


P 2 

a + g + Y 

n,m  m,n  m,n 


= 1 


^m.n"  A _ Sin  <}>)  - 7 (y  Cos  <J>) 

m m ,n  m 


+k(g_  Cos  <*>— ot  Sin  <J>)} 

m>n  m 


(11,3) 


Similarly  to  express  ^m.n  ^m.n  x Tm,n  we  proceed  in  the  above 
manner  to  get 

nm,n=  ~ ([(yf  _ + 3^  ) Cos  <p-a  B Sin<f>]  i 

n m,n  m,n  m,n  m,n 


+ [(y!I  r,  + at  ni  Sin,fr  ~am  Cos*>  ] J 

m,n  m,n J m,n  m,n 


-(y_  _ & _ Sin <J>  +y  _ a Cos*]  k} 
m,n  m,n  m,n  m,n 

4.  Hydrodynamic  Coordinates  for  the  Tie  Leg 

The  tie  leg  is  a single  cable  and  as  such,  does  not  involve  an  index 

n,  with  the  results  the  expressions  (10,3),  (11,3)  and  (12,3)  are  valid  if 
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Thus  the  transformation  expressions  for  hydrodynamic  coordinates  in  case 


of  a tie  leg  would  he  expressed  as 


5.  The  Hydrodynamic  Distributed  Forces 


Once  a set  of  hydrodynamic  coordinates  has  been  obtained,  it  becomes 


a simple  matter  to  represent  the  hydrodynamic  forces  for  unit  length  which 


are  then  resolved  into  x,  ir  and  n directions 


inherent  difficulties  that  exist  in 


defining  hydrodynamic  force 


Limited  data  is  available  to  determine  tne 


functional  forms  of  a side  force  in  the  tt  direction.  However,  it  is  known 


that  magnitude  of  this  force  compared  with  the  other  two  components  is  small 


and  as  such  is  neglected 


T 

m,n 


A 

m,n 


A is  the  sine  of  the  angle  between  the  (m,n)th  cable  segment  and  the 
ID  y n 


stream  and  is  given  by 
A 


m,n-*^  2 

y + 


Y_  _ + (3  Cos  *-a  Sin  <|>  )‘= 
m,n  m,n  m,n 


so 


I = V(Z) 


Tm,n  J~ 2 

n 


2 2 

a + 3 + Y 

n,m  m,n  'm,n 


= 1 


as 


V.n-  A {1W  „ Sin  <(>)  - J (y  Cos  *) 

m,n  m,n  m,n 


+k(B  Cos  *-a  Sin  *)} 
m,n  m,n 


(11,3) 


4.  Similarly  to  express  ^mjn  =1Tm,n  x Tm,n  we  proceed  in  the  above 
manner  to  get 


'V.n3  {[(y2  + g2  ) Cos  4>— oc  (5  Sin<f> ] i 

* A m,n  m,n  T m,n  m,n 


+ „ + rO  Sin*  ~am  nPTn  nCOS()>  ] J 

m,n  m,n;  m,n  m,n 


-t Y 3 Sin*  +y  o Cos*]  k} 
m,n  m,n  m,n  m,n 


Hydrodynamic  Coordinates  for  the  Tie  Leg 


The  tie  leg  is  a single  cable  and  as  such,  does  not  involve  an  index 
with  the  results  the  expressions  (10,3),  (11,3)  and  (12,3)  are  valid  if 


b)  Normal  Drag  Force  in  n Directions 


It  is  known  that  the  normal  drag  force  for  unit  length  which  acts  in 


p = the  density  of  the  fluid 

= the  coefficient  of  drag  of  (m,n)th  cable  segment 
i,n 

= the  diameter  of  (m,n)th  cable  segment 


Expression  for  normal  drag  force  in  n direction  for  the  tie  leg  array 
is  essentially  the  same  as  (l6,3a)  except  that  no  index  n exists  and  as  such 


In  most  of  the  work  done,  this  component  is  neglected  being  made  equal 


this  analysis  assume 


that  this  force  can  be  given  by  the  expression 


the  coefficient  of  drag  of  (m,n)th  cable  segment  when  the  segment 


where  (^m  is  coefficient  of  drag  of  mth  cable  segment  when  the  segment  is 
parallel  to  the  stream. 

6.  The  Total  Hydrodynamic  Force 

From  the  above  expressions,  then,  it  becomes  clear  that  the  hydrodynamic 
force  per  unit  length  which  is  acting  on  the  (m,n)th  cable  segment  can  be 
written  as: 


fm,n  fm,n  nm,n  fm,n  xm,n 


(18,3) 


where 


n 

fm,n 


= - CN  d (V(Z)A  y 

2 m,n  m,n  m,n 

= yC  .2  V2(Z) 

m,n  Am,n 


(19,3) 


and 


yC  -fc*  d 
m,n  2 m,n  m,n 


(20,3) 


also 


f'  4c-  d V(Z)  [V.Tm,n] 
m,n  e m,n  m,n 


=yf  _ r^  V2(Z)[a  Cos  <(>+$  Sin$] 


m,n  m ,n 


m,n 


m,n 


(21,3) 


where 


D _ CT2 m, 

" c’». 


(22,3) 


Thus  having  known  the  expressions  for  fn  m.n  and  fT  m,n  - these  can  be 
substituted  into  the  equation  (l8,3)  to  get  the  total  hydrodynamic  force. 
Similarly,  the  expression  for  total  hydrodynamic  force  in  the  case  of.  a tie  1«g 
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is  given  by 


Total  n . _t 
f = f n + f t 

m m m ram 


(23,3) 


where  f*  and  fn  are  similar  to  Eqs.  (21,3)  and  (19,3)  with  no  (n) , 
in  m 

To  derive  expressions  for  the  projection  of  these  hydrodynamic  force 
in  the  i,  J and  k dirn  , we  proceed  as  follows: 

The  hydrodynamic  force  for  the  tri-moored  structure  acting  in  the  a dir11 
is  given  by 


i 

i 


Cc.x  v*(z) 
m,n 


(23,3a) 


Similarly 


hc;y  = Cc:y  jt{) 
m,n  m,n 


(23,3b) 


hc:z  = Cc:z  V^Z) 
m,n  m,n 


(23,3c) 


and  for  the  tie  leg  by 

hC:X  = c"*  V^Z) 

m m 


(24,3a) 


hC:y  = cc:y 


v2(z) 


(24,3b) 


hC:y  = CC:y  V^Z) 


(24,3c) 


where  C m,n's  are  known  as  the  drag  constants  and  they  are  defined  by: 

„c:x  c r.  / 2 „2  v „ 

C _ = u [A  (y  + 6 ) Cos  <t>  - 

m,n  m,n  m,n  m,n  m,n 


1 


A g Sin  + 

m,n  am,n  m,n 


rm,n  (“m,n  Cos*  + 6m,n  Sin*}  ° m’n]  (?5’3a) 


■ 


and 


Cc:y  = pC  [A  (y2  + ot2  ) Sin*- 

m,n  m,n  m,n  rm,n  m,n 


A o 6 Cos*  + 
m,n  m,n  m,n 


rm  r»  (“n,  n COB<fr  + n Sin<fr)&„  „ 1 

m,n  m,n  m,n  m ,n  J 


(25,3b) 


Cm  n = „[-Am  „ (“  „ Cos*  + B Sin<|> ) Ym,n 

m,n  m,n  m,n  m,n  m,n 


+ „(“m  „ Cos*  + B „ Sin*)Ym,n] 

m,n  m,n  m,n 


(25,3c) 


similarly  for  the  tie  leg 


C™  X = + ef)  Cos  * -A  a 8 Sin  *+ 

m m m m m mmm 


r (a  Cos*  +8  Sin*)am] 
mm  m 


(26,3a) 


''m  y = y*(An,  (y^  + af)  Sin  * -A  a 8m  Cos  * + 
m mmm  m mmm 


r»  (“-,  Cos(t'  +B  Sin*) 6 ] (26,3b) 

m m m m 

and 

C * Z C r D 

C = u [ -A  (a  Cos*+8  Sin*)Y  + r 
m mmm  m mm 

(a  Cos*  + 8_  Sin  *)y_]  (26,3c) 

m m in 

Once  the  hydrodynamic  force  on  the  cable  array  is  found, the  next  stage 
is  to  lump  these  forces  on  the  stations.  To  lump  the  distributed  forces  as 
a single  force  which  acts  at  the  (m,n)th  cable  station,  the  half  segment 
lumping  technique  (5)  is  employed.  Then  using  this  technique,  the  lumped 
hydrodynamic  forces  are  given  b>  the  equations: 


UO 


for  m = 2 ,3 M(n)-1  n = 1,2,3 


Ed0  - CC:6  V (m,n:  BL  /0,  BL  ) + 
m,n  m,n  m,n/2’  m,n 

£l,n  V (nH,n:  0,  ^^1 

and  for  m * M(l)  and  n = 1 


'SfiD.l  “ °M(l),l  7 BLH(1),1/2,BLM(1),1^ 

+ °J;l7<1>2i  °-BL1.2/2> 

♦ cj;e3  7(1,2:  0,  bl1i3/2) 


,c:0  77 


(27,3) 


(28,3) 


where  6 = x,y  or  z 

The  function  V (m,n,|  jj)  represents  the  integral  of  (V  (Z)  )along  the 
(m,n)th  half  segment  and  is  represented  by 

K 

(29,3' 


V (m,n,5  ,5  ) = (v2[Z(m,n :£ ) ]d£ 
± * 51 


where  the  argument  z of  V (Z)  is  expressed  in  terms  of  the  integration 
parameter  5 along  the  (m,n)th  cable  segment  through  the  relation 


Z = Z (m,n:£)  = Z + y 

m,n  m+l,n 


(30,3) 


where  ym+1  q is  the  direction  cosine  of  the  mth  segment  as  defined  oy 
equation  (10,3c). 

In  the  case  of  a tie  leg  the  half  segment  technique  yields  the  lumped 
hydrodynamic  forces  given  by  the  following  expression: 


1*1 


for  m = 2,3 MN-1 

Hc:6  = Cc:0  V (m:B  . ,B  )+ 
m m m/2  ’ m 


m 

,c:6 


-»*i  v t"*1- 0 • Vi>  <31'3> 

vbere  0 = x,  y or  z 

and  expression  V (m,El,S2)  represents  the  integral  of  V^(z)  along  the  mth 

half  segment  and  is  represented  by 

V (m,S  5 ) = Pv2  [Z(m,£) ] d£  (32,3) 

i,  2 / 

"1 

D.  DISCRETE  ELEMENT  HYDRODYNAMIC  FORCES 

The  direction  of  drag  on  the  discrete  elements  becomes  parallel  to  V if 
it  is  assumed  that  there  is  no  lift  associated  with  them.  The  magnitude  of 
the  drag  force  in  that  case  is  given  by: 

P 


- C?  A V^Z) 

2 k,m,n  x,m,n 


(33,3) 


where 


C^  m n “ coefficient  of  drag  of  the  (k,m,n)th  elemental  device 
A,  = effective  cross-sectional  area  of  (k,m,n)th  elemental  device 
Then  if  this  force  be  resolved  in  x,  y and  z coordinate  system,  the 
hydrodynamic  force  due  to  the  (k,m,n)th  elemental  device  is  found  as: 

C.n  * C.u  »*»>  <3M*> 


h = c.e:y 

K.9m9n 


X ,m,n 


v^z) 


(34,3b) 


e:z 


^.n  = 0 (34,3c) 

where  the  drag  constants  of  the  (k,m,n)th  elemental  device  are  given  by 

(35,3a) 


C*:X  = v?  Cos$ 
k,m,n  k,m,n 


„e  :y  e , 

k,m,n  Mk,m,n  in* 

0 

and  y,  being  the  hydrodynamic  coordinate  is  defined  by 
&.  ,m  g n 

e s t c D a 

^k,m,n  2 k,m,nx,m,n 


(35,3b) 


(36,3) 


Similarly,  for  the  tie  leg  the  hydrodynamic  force  due  to  the  (j,m)th  elemental 


device  is  given  as 

he;x  = ce:x  y2(z) 
j,m  J,m 


(37,3a) 


he:y  _ ce:y  y2(z) 
J »® 


(37,3b) 


h1m  = 0 

j ,m 


(37,3c) 


where  drag  constants  of  the  (j,m)th  elemental  device  are  given  by 


_e:x  e 
CJ,m  " ^ Cos* 


(38,3a) 


C1:I  " A m Sin* 

J ,m  J ,m 

0 

and  y is  the  hydrodynamic  coordinate  defined  by 
J 


(38,3b) 


e = P D 
u j ,m  2 j ,m  m ,n 


(39,3) 


The  next  step  is  to  apply  the  half  segment  lumping  technique  to  find  the 
hydrodynamic  forces  acting  at  the  (m,n)th  station  and  due  to  discrete 
elemental  devices.  This  is  given  through  the  following  relations: 
for  m = 2,3  ...M(n)  - 1 and  n = 1,2,3 


iie:6  _ k(m,n)  E:0  ,_2  r_,  _ .. 

v C.  » LZ(m,n :S,  )1 

m,n  k=I(m,n)+l  k»m»n  k»m’n 


k(m+l,n) 

, C,e* 


k = 1 k,uH-l,n 
where  0 = x,  y or  z 
and  for  m = M(l)  and  n = 1 


Mi.  v^[Z(m+l,n:S.  - )] 


k,m+l,n 


( 1*0 ,3) 


1*3 


e:0  _ k(M(l),l) 

»(l).l  " JL  CkM(l)  iv  [z(M(l).1:f  m(D  JJ 

k=k(M(l)  ,l)+l  k’MU,»1  E>bU,>1 


k(l»2)  e-0  2 

*(l,3)  .-8  2 

c;;!,3  ^ tz(i‘3:sk>i.3)1  (U1'3) 

Again,  6 = x,  y or  z and  argument  z of  (V2(Z))  is  expressed  in  terms  of 
the  position  of  the  (k,m,n),  the  device  along  the  (m,n)th  segment  through  the 


relation 


Z(m,n:  S.  ) = Z , + y S. 

k,m,n  m-l,n  m,n  k,m,n 


(42,3) 


Using  half  segment  technique  to  find  hydrodynamic  forces  acting  at  the 
mth  station  of  the  tie  leg  due  to  discrete  elements,  the  following  expression 
is  derived: 

II®'®  — j(®)  » • A 2 

m " I C*-®  [V2  Z(m,S  )] 

J=y(m)+1  J »m  J *m 


+ J(n+1)  e-8  2 

£ C i'Li  ^ S,  - ) ) 

is  i J*®*!  J ,m+l 


(43,3) 


for  m + 2,3  MN-1 

where  9 * x,  y or  z 


E.  FINAL  EXTERNAL  FORCES 


Thus,  to  summarize  the  external  force 

(C- 

that  is  acting  at  the  (m,n)th  cable  station  for  the  tri -moored array  can  be  given 


as : 


and 


F*  = 

Hc;x 

+ 

HejX 

m,n 

m,n 

m,n 

& * 

Hc;y 

+ 

H® ;y 

m,n 

m,n 

m,n 

FZ 

W 

+ 

HejZ 

m,n 

m,n 

m,n 

(1*1+ ,3a) 
(UU,3b) 


(1+1+ ,3c) 


Where  the  lumped  weight  forces  W are  defined  by  Eqs.  (5,3),  the  lumped 

xu  ,n 

c * 0 

hydrodynamic  forces  H * due  to  drag  forces  on  the  cable  segments  are  defined  by 

m ,n 

e * 0 

Eqs.  (27,3),  and  the  lumped  hydrodynamic  forces  H ’ due  to  the  drag  forces  on 

m ,n 

the  elemental  devices  are  defined  by  Eqs.  (1+0  ,3) . The  above  equations,  however, 
do  not  apply  to  stations  (M(2)  ,2)  and  (M(3),3)  as  imaginary  reactions  act  on  these. 
Similarly,  the  external  force 


(F*.  J*.  FZ) 
m*  m*  m 


acting  at  the  m th  cable  station  for  the  tie  leg  array  can  be  given  by  the  following 
equations : 

(U5,3a) 


( 1+5  ,3b) 


1^- 

HC;*  + 

He  ’X 

m 

m 

m 

**  = 

HC;y  + 

He;X 

m 

m 

m 

and 

F^  = W + H°;z  (1*5, 3c) 

m m m 

Where  W^  is  defined  by  Eqs.  (l+,3),  is  defined  by  Eqs.  (31,3)  and 

H6’0  is  defined  by  Eqs.  (1+3,3). 
m 

In  the  above  equations,  only  the  weight  forces  are  constant.  The  hydro- 
dynamic  forces  depend  on  the  position  of  the  structure  through  both  the  orienta- 
tion and  the  depth.  As  discussed  in  Chapter  I,  this  problem  is  solved  by  using 
the  Method  of  Imaginary  Reactions  in  conjunction  with  the  Method  of  Successive 
Approximations . 


1+5 


CHAPTER  IV 


THE  CONVERGENCE  ALGORITHM  USING  A BINARY  SEARCH  ROUTINE 

In  Chapter  II  it  was  suggested  that  a convergence  algorithm  is 
required  to  satisfy  the  force  balance  condition  using  a near-minimum 
number  of  subroutine  iterations.  The  force  balance  constraints  suggest 
that  in  order  to  obtain  the  tie  leg  effect  in  a lambda  structure  (refer 
to  Fig.  6)  the  following  equations  must  be  satisfied. 

D = D <1'U) 

E ’ ral*.t  E 

where  FFY  and  FFY*  are  the  resultant  forces  exerted  by  the  main  cable 
arrays  on  the  tie  leg  and  RY1  and  RY1*  are  the  end  reactions  that  the 
tie  leg  exerts  on  the  main  cable  array. 

If  the  above  set  of  forces  do  not  satisfy  the  force  balance  constraints , 
then  a different  value  for  FFY  and  FFY*  is  used  which  produces  a different  chordal 
distance  D’E’  such  that  the  tie  leg  fits  into  this  new  distance.  Now, 
the  force  balance  test  is  applied  and  if  it  is  not  satisfied,  the  above 
process  is  repeated. 

This  chapter  discusses  a technique  which  determines  a new  value 
for  forces  FFY  and  FFY*  so  that  the  number  of  iterations , before  forces 
converge  to  satisfy  equations  (1,U),  is  hopefully  made  minimal.  The 
following  discussion  is  restricted  to  force  FFY  and  reaction  RY1;  however, 
the  same  discussion  holds  for  force  FFY*  and  reaction  RY1*. 

Various  algorithms  to  find  a new  value  of  FFY  (refer  to  Figure 
6)  were  evaluated;  however,  a Binary  Search  Routine  was  found  to  be 
the  most  direct  method  to  achieve  fact  convergence. 
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A.  BINARY  SEARCH  ROUTINE 


A Binary  Search  Routine  is  an  algorithm  used  to  find  a new  value 
of  FFY  such  that  it  always  lies  in  between  the  previous  two  values  of 
FFY  and  RY1.  If  the  force  FFY  produces  a reaction  RY1  in  the  tie 
leg,  then  the  new  value  of  FFY  is  given  by: 


FFY 1 =ABsj  ABS1  RYIj[-ABSiFF-Yi|  + Minimum  (ABs|fFy|  ,ABs|RYl|)  (2,4) 

This  new  value  of  FFY  as  given  by  FFY'  and  used  with  the  proper  sign 
on  the  tri-moored  array,  then,  produces  another  value  of  RY1.  The 
force  balance  criterion  is  applied  and  the  process  repeated  until  the 
condition  is  satisfied.  A typical  plot  of  the  result  of  using  this 
algorithm  in  a practical  problem  is  as  shown  in  Figure  10,  where  the 
convergence  can  be  seen. 


— ■■■  --  - — . noj  .of  iterations 

Figure  10:  The  Behavior  of  the  Force  FFY  When  Binary  Search  Routine 

Is  Used 

The  details  of  the  Binary  Search  Routine  as  used  in  the  problem  are 


listed  below. 


If  FFY  is  the  first  assumed  force  that  is  applied  on  the  main  cable, 
then  after  equilibrium  has  been  attained,  a chordal  length  into  which  the 
tie  leg  is  to  fit  is  obtained.  If  this  equilibrium  configuration  of  the 
main  structure  is  not  compatible  for  the  tie  leg,  then  another  assumed 
value  FFY ' (FFY'  > FFY)  is  used  on  the  main  cable  which  reduces  the 
chordal  distance  by  an  amount  X ( see  Figure  11) . 

Thus  by  several  iterations  a tie  leg  is  made  to  fit  into  an 
appropriately  adjusted  chordal  distance,  as  shown  in  Figure  12.  This 
produces  a reaction  RY1  in  the  tie  leg.  Now  the  force  balance  test 
is  applied  and  if  it  is  not  satisfied,  a new  value  of  FFY" (FFY'<  FFY"<  RYl) 
is  obtained  from  the  Binary  Search  Routine  using  equation  (2,4).  This 
process  is  repeated  until  the  force  balance  condition  is  met.  When 
this  happens  equations  (1,4)  have  been  satisfied. 


The  force  FFY , drawn  as  a function  of  x is : 


Similarly  if  a curve  of  reaction  RY1  is  drawn  with  respect  to  X,  it  will 


If  both  these  curves  are  superimposed,  a graph  as  shown  in  Figure  13  is 
obtained. 

Then  to  summarize,  FFY'  produces  the  reaction  RY1  and  this  is  represented 
by  AB  (notation  follows  Figure  13).  Using  binary  search  a new  value  of  FFY 
is  found  and  is  represented  by  CD.  Force  FFY  at  D produces  a reaction  RY1 
given  by  E so  that  this  case  is  represented  by  DE.  Binary  Search  Routine  is 
used  again  to  find  the  new  value  of  FFY  represented  by  FL. 

Force  represented  by  FL  is  used  to  find  the  chordal  distance  in  which 
the  tie  leg  has  to  fit.  A reaction  represented  by  point  Pis  required  to  do 
so.  The  Binary  Search  Routine  is  again  used  to  find  a new  value  for  force 
FFY,  which  is  now  represented  by  MN.  Thus,  this  process  is  repeated  until 
the  forces  represented  by  point  K are  obtained.  At  this  stage  equations  (1,4) 
are  satisfied.  The  convergence  to  the  required  forces  that  satisfy  the  force 
balance  condition,  in  all  computer  runs  and  tests  made,  appears  to  be  very 
qui ck . 

A block  diagram  representing  the  stepwise  use  of  the  Binary  Search 
Routine  is  shown  in  Figure  1^ . 


FORCES 


Apply  the  forces 
(FFX,FFY & FFZ) 
on  the  main  array 
to  the  tie  leg 


This  will  result  in  a chordal 
distance  into  which  the  tie 
leg  has  to  fit 


Subroutine  tie  leg 
will  yield  forces 
(RX1 , RY1  & RZ1 ) 
used  to  fit  the 
tie  leg  into  the 
chordal  distance 


Apply  Force  Balance 
criteria  to  FFY, 

& RX1 


IF  IT  IS  SATISFIED 


IF  IT  IS  NOT  SATISFIED 


Use  Binary  Search 
Routine  to  find 
another  value  for 
FFX,  FFY,  FFZ 


OUTPUT 


Figure  : Block  Diagram  of  the  Binary  Search  Routine 
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B.  PRECISION  FOCUS 


As  described  in  Chapter  I , the  equilibrium  configuration  of  the  array- 
system  is  obtained  if  the  value  of  E,  the  measure  of  error,  is  nearly  zero. 
Theoretically,  the  iteration  could  continue  until  E is  exactly  zero.  However, 
this  is  unnecessary  for  useful  accuracy  and  a cut  off  value  COMPE  is  defined 
that  determines  the  acceptable  completion  of  the  iterative  process.*  In  the 
course  of  developing  the  computer  analysis  it  was  found  that  the  tie  leg  array 
would  not  converge  to  its  equilibrium  configuration  within  specified  limits 
of  COMPE.  This  was  attributed  to  the  higher  levels  of  required  forces  in 
the  y-coordinate  direction  as  compared  to  the  other  two  directions,  with  the 
result  that  6 becomes  too  small  for  rapid  convergence  in  all  three  directions. 


To  overcome  this  problem,  another  convergence  technique,  called  Precision 

Focus,  was  used  in  conjunction  with  the  method  of  convergence  that  has  been 

(13) 

described  in  Chapter  I.  Precision  Focus  has  been  used  by  Savage  and  is 


.1 
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CHAPTER  V 

TESTING  THE  COMPUTER  MODEL 

The  computer  program  written  to  simulate  the  given  array  system  is 
listed  in  Appendix  I.  This  chapter  presents  a discussion  of  some  testing 
of  this  computer  model,  and  some  aspects  of  the  general  behavior  of  the  array 
system. 

The  computer  model  tested  has  the  following  specifications: 

The  tri-moored  structure  consisted  of  three  identical  legs, 
each  25,000  ft.  long,  in  the  unstressed  state.  These  legs  have  a 
diameter  of  0.675  inches  and  a weight  in  water  of  O.606  lh/ft. 

The  extensional  rigidity  of  each  cable  segment  is  given  as  2.2x10^ 
lb.  The  co-efficient  of  normal  drag  is  assumed  as  1.U0  for  the  entire 
range  of  current  velocities. 

The  tie  leg  cable  has  a diameter  of  0.675  inches  and  an  extensional 
rigidity  of  2.2x10^  lb.  The  length  and  weight  of  the  cable  was 
determined  internally  by  the  program. 

As  described  in  Chapter  II,  the  basic  approach  in  the  analysis  of  the 
given  system  has  been  to  divide  the  array  system  into  the  following  coupled 
parts : 

1.  The  tri-moored  array  structure 

2.  The  tie  leg  array  structure 

This  approach  was  thought  to  be  most  practical  because  results  obtained,  namely 

the  displacements  of  the  subsurface  buoy  in  the  given  system, could  be  compared 

readily  with  the  displacements  obtained  from  the  analysis  of  a tri-moored  structure 

(9) 

without  a tie  leg  for  which  a computer  program  has  been  published.  Also, 
it  is  possible  to  study  the  behavior  of  the  main  structure  and  the  tie  leg 
separately.  To  check  the  program  output  the  following  extreme  cases 
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were  analyzed  on  the  computer. 

1.  When  the  tie  leg  is  attached  near  the  top  of  the  tri-moored  structure, 
where  the  tie  leg  has  negligible  length 

2.  When  the  tie  leg  is  attached  near  the  bottom  of  the  tri-moored  structure, 
where  the  tie  leg  is  essentially  connected  between  two  of  three  anchor 
points 

Since  both  of  these  are  limiting  conditions  where  the  tensions  will  be 
taken  up  by  the  anchors  in  one  case  and  where  the  tie  leg  is  of  negligible 
length  in  the  other , the  system  should  behave  similarly  to  a simple  tri-moored 
structure  as  shown  in  Figure  1.  To  look  at  these  extreme  conditions  and  at 
the  general  case  of  the  tie  leg  attached  anywhere  on  the  main  structure,  the 
computer  model  was  used  and  the  results  are  discussed  in  the  next  three 
sections. 

A.  TESTING  THE  COMPUTER  MODEL  WHEN  THE  TIE  LEG  IS  NEAR  THE  TOP  OF  THE  TRI-MOORED 
STRUCTURE 

The  computer  program  is  written  in  such  a manner  that  it  is  not  possible  to 

place  the  tie  leg  precisely  at  the  apex.  To  cope  with  this  problem,  it  was  decided 

to  break  the  main  cable  into  unequal  segments  such  that  the  first  and  last  segments 

were  only  of  5 feet  length.  Each  of  the  cables  was  broken  up  into  20  segments - 

This  resulted  in  a simulation  which  allowed  the  tie  leg  to  be  put  at  a distance 

of  five  feet  from  the  top.  The  displacements  obtained  from  this  configuration 

were  used  for  comparison  purposes.  Many  different  conditions  of  tie  leg  parameters 

were  tried  and  in  each  case,  the  tie  leg  structural  analysis  gave  similar  results 

(g) 

to  the  already  tested  program  for  a simple  tri-moored  structure  . For  example, 
with  a tie  leg  of  6.12  feet  in  length  and  the  weight  per  foot  of  cable  at  0.505 
lb. /ft.,  the  displacement  of  the  buoy  due  to  a standard  test  current,  for  the 
case  of  -arTri-moored  array  without  this  tie  leg  is : 

51* 


Horizontal  deflection  of  the  buoy  = 36.95  ft. 


f 
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Vertical  deflection  of  the  buoy  = -11.50  ft. 

The  displacement  of  the  buoy  due  to  the  same  test  current  with  the  tie  leg 
cable  attached  between  two  cables  near  the  apex  of  the  tri-moored  system  is : 
Horizontal  deflection  of  the  buoy  = 35.12  ft. 

Vertical  deflection  of  the  buoy  = 11.50  ft. 

The  results  obtained  from  this  test  condition  compare  favorably  with  the  results 
obtained  from  the  analysis  of  a simple  tri-moored  structure. 

B.  TESTING  THE  COMPUTER  MODEL  WHEN  THE  TIE  LEG  IS  LOCATED  NEAR  THE  BOTTOM  OF 
THE  TRI-MOORED  STRUCTURE 

For  this  condition,  the  tie  leg  was  attached  near  the  bottom  (five  feet  up 
each  leg),  between  two  of  three  anchor  points. 

The  length  of  the  tie  leg  was  30612.8  feet  and  weight  of  the  cable  per 
foot  was  0.367  lb.  This  weight  of  the  cable  was  determined  internally  by  the 
program  so  as  to  make  the  tie  leg  array  slightly  positively  buoyant. 

The  displacements  of  the  buoy  when  placed  in  the  test  current  with  the 
tie  leg  are  as  follows : 

1.  Horizontal  Deflection  of  the  buoy  = 39*61  ft. 

2.  Vertical  Deflection  of  the  buoy  = -12.1+5  ft. 

These  results  should  be  compared  with  the  displacements  of  the  standard 
tri-moored  buoy,  as  listed  above. 

C.  PRELIMINARY  STUDY  USING  THE  COMPUTER  MODEL 

Besides  the  above  two  test  cases , the  installed  position  of  a tie  leg  was 
varied  along  the  length  of  the  main  cables.  This  was  done  in  order  to 
study  the  behavior  of  the  model  for  different  positions  of  a tie  leg. 

While  the  purpose  of  this  report  is  not  to  conduct  a complete  study  of 


system  'behavior,  enough  data  was  gathered  to  predict  the  general  behavior 
of  a tie  leg  in  a tri-moored  structure.  This  behavior  is  presented  in 
Figure  15  and  l6  wherein  the  vertical  horizontal  deflections  of  the  buoy 
are  plotted  against  the  position  of  the  tie  leg.  The  curves  should  be 
considered  only  as  approximate  since  the  data  points  were  obtained  from 
calculations  using  slightly  different  tie  leg  cable  weights  for  each 
location.  This  was  done  as  a matter  of  convenience.  The  error  introduced 
does  not  alter  the  general,  behavior  or  the  conclusion  that  horizontal 
deflections  of  the  apex  are  increased  by  approximately  a factor  of  two 
and  the  vertical  deflections  of  the  apex  are  increased  by  approximately 
a factor  of  four.  If  a prototype  system  is  contemplated,  a more  detailed 
analysis  should  be  conducted. 

From  the  graphs  it  can  be  inferred  that  worse  deflections , both 
horizontal  and  vertical,  are  encountered,  as  one  would  expect,  when  the 
tie  leg  is  at  a position  near  the  middle.  The  current  used  is  perpendicular 
to  the  tie  leg. 
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This  curve  should  be  • \ 
considered  only  as  spprsriwste  \ 
since  the  data  paints  wars 
obtained  flww  calculations 
using  slightly  different  tie 
leg  cable  weights  far  each  location 


Position  of  tie  leg  along  the  main  cable  (station) 

>:  The  Effect  of  the  Position  of  the  Tie  Leg  on  Horizontal 


o 

Deflection  of  the  Buoy 
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This  curve  should  be  considered 
only  as  approximate  since  the  data 
points  were  obtained  from  calculations 
using  sll^itly  different  tie  leg 
cable  weight  for  each  location. 


0 2 3 4 5 6 7 8 9 10  tl  12  13  14  15  16  17  18  19  20  2) 

Position  of  the  tie  leg  along  the  main  cable  (station) 

Figure  16:  'Hie  Effect  of  the  Position  of  the  He  Leg  on  Vertical 
Deflections  of  the  Array 


APPENDIX  I 


THE  COMPUTER  PROGRAM 


“ 


A.  DESCRIPTION 

The  computer  program,  required  to  simulate  a tri-moored  buoy  with  a tie 
leg  is  reproduced  in  this  appendix.  This  program  is  written  in  FORTRAN  IV 
and  can  be  compiled  and  executed  on  most  of  IBM-360  facilities.  The  facility 
at  the  University  of  New  Hampshire  is  IBM  360-1+0  and  the  computer  program 
has  been  specifically  adapted  to  it. 

The  program  consists  of  lU  subroutines  besides  the  main  section.  Comment 
cards  at  the  beginning  of  each  of  these  sections  describe  the  nomenclature 
and  in  some  cases  the  purpose  of  each  subroutine. 

The  program,  as  written,  is  restricted  to  20  segments  and  to  10  elemental 
devices  per  segment  in  each  of  the  main  arrays  and  to  21  segments  and  5 elemental 
devices  in  each  . tie  leg  array.  This  is  only  for  convenience  and  these 
numbers  can  easily  be  changed  by  changing  the  dimensions  of  the  common  arrays 
in  the  main  program  and  in  the  subroutine  and  functions. 

Since  the  length  of  the  tie  leg  is  dependent  upon  the  pretensioning,  the  length 
of  each  segment  is  determined  internally.  The  tie  leg  is  made  positively 
buoyant  by  about  2%  and  as  such,  the  weight  per  unit  feet  is  also  computed 
internally. 

The  total  length  of  the  computer  program  is  (FBJjS)^.  If  at  any  facility 
the  computer  memory  is  inadequate,  then  K and  M dimensions  for  all  the  arrays 
are  reduced  to  the  largest  values  Of  KMAS  (m,n)  and  JMAX(m)  in  the  case  of  the  tie 
leg — depending  upon  any  particular  analysis. 

B.  INPUT  DATA  CARDS 

The  complete  input  data  is  controlled  by  subroutine  INPUT.  Thus,  data  cards 
in  the  program  correspond  to  this  subroutine. 
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This  subroutine  has  been  broken  up  into  four  parts: 

First  part  consists  of  reading  in  of  data  that  is  valid  for  the  main 
arrays  and  the  tie  leg  arrays. 

Second  part  deals  with  data  relevant  only  to  the  main  arrays. 

Third  part  consists  of  reading  in  of  data  for  the-  tie  leg  array. 

Finally,  the  last  part  deals  with  the -profiles  of  v.Le  hydroayrun.lc  current  as  n 
attacks  the  main  array  and  the  tie  leg  respectively. 

1.  The  first  input  card  contains:  COMPE,  COMPD,  STAPSI , DELPSI , ENDPSI,  and 
TIECOM.  F10.3  FORMAT 


COMPE 

COMPD 

STAPSI 

DELPSI 

ENDPSI 

TIE  COM 

Comparison 
value  for  E- 
Error  Function 

Comparison 
value  for 
displace- 
ment 

First  cur- 
rent angle 
to  be 
i analyzed 
in  leg 

Change  of 
current 
angle  in 
Idegrees 

Final 

current 
angle  in 
r degrees 

Comparison 
value  to 
meet  force 
balance 
criteria 

2.  The  next  three  cards  contain  the  Anchor  Positions,  AA1  = (=X)  , BB1  (=Y), 
CC1  (=Z)  -F  10.3  FORMAT 


card  2 
card  3 
card  b 


3.  The  fifth  card  contains  the  number  of  stations  per  cable,  MMAX(N)  - I 5 FORMAT 


AAl(=X) 

BB1(=Y) 

CC1(Z) 

X of  Anchor  1 

X of  Anchor  2 

X of  Anchor  3 

1 

Y of  Anchor  1 

Y of  Anchor  2 

Z of  Anchor  3 

Z of  Anchor  1 

Z of  Anchor  2 

Z of  Anchor  3 

MMAX(l) 

MMAX(2) 

MMAX(3) 

No.  of  stations 
on  cable  1 

No.  of  stations 
on  cable  2 

No.  of  stations 
on  cable  3 

b.  The  next  group  of  cards  contain  the  physical  properties  of  the  cable  segments 
and  the  discrete  elements.  This  starts  with  cable  1,  which  goes  from  anchor 
point  to  subsurface  buoy — followed  by  cables  2 and  3. 


FORMAT  OF  CARDS  DESCRIBING  THE  PROPERTIES  OF  THE  CABLE  SEGMENTS 


BLBAR 

wc(m,n  ) 

xten(m,n) 

tdrag(m.n) 

(M,N) 

F 10.2 

F 10.2 

F 10.2 

F 10.2 

Unstressed 

length 

(ftj 

Weight/ 

ft. 

(lb. /ft.) 

Extensional 

regidity 

(lb-) 

Normal  drag 
coefficient 

CABDIA 

(M,N) 

F 10.2 

PDRAG(M.N) 

F 10.2 

KMAX(M,N) 

1-5 

Cable 

Dia. 

(inches ) 

Parallel 

Drag 

Coeff. 

No.  of 

Discrete 

Elements 

following  each  one  of  these  is  KMAX(m,n)  cards  giving  the  physical  properties 
of  the  discrete  elements 


SBAR(k,m,n) 

WE 

(k,m,n) 

DRAG  CF(k,m,n) 

XREA(k,m,n) 

Length 

weight 

Drag,  coeff. 

X-sectional 

segment 

lb. 

area 

(m-1 ,n) 

element 

5.  The  next  group  of  data  cards  contain  input  data  for  the  tie  leg  array. 
1.  First  card  in  this  group  contains  the  number  of  station  on  the  main 
array  to  which  tie  leg  is  attached,  - I 3 FORMAT 


1-13 

No.  of  stations  on  the 
main  array  to  which  tie 
leg  is  attached 


'ta 


1 


2.  The  next  consists  of  number  of  stations  into  which  tie  leg  is  broken 
I 2 FORMAT 


MMAXN 

1-2 


No.  of  stations  on 
tie  leg 


3.  The  next  group  consists  of  data  for  the  physical  properties  of  the 
tie  leg  array  segments. 


XXTENiM 


Extensional 

rigidity 

(lb.) 


TTDRAGIM 


Normal 
drag 
coeff . 


cabda(m 


Cable 
diameter 
( inches ) 


pfdrag(m) 

JMAX(M) 

Parallel 

drag 

coeff. 

No.  of 

discrete 

elements 

following  each  one  of  these  is  JMAX(M)  cards  giving  the  physical  properti 
of  the  discrete  elements 


WEE(J,M) 


DRAGCF(J.M) 


Weight  (lb.)  Drag  coeff. 


XXREA(J,M) 

X-sectional 

2 

area  ft. 


Finally  the  last  group  of  data  cards  deal  with  the  input  values  for  the 
current  profile. 


H(k)  F 10.3 

V(k)  F 10.3 

Z coordinate  where  the 
velocity  profile  changes 
slope  (ft.) 

The  magnitude  of 
current  at  Z=H(k) 
(fVsec. ) 

nor*  f* o non  i : ^n*  i ofio 


COMMON/ Cl 9/ T I ECOM,  L 


GALL  INPUT 


ikicnou  . firm 

■ If*  * 


GGHPUTE  *IOSEGEHENT  DI-SGRETE  ELEMENT 


MX  = MMAX(N) 

XX  * KMAX(MyN) 

I F( S BAR! K»M»N I . GT .BL BAR (M«N)/2. IGO  TO  2 

i rntir  tkinc 


K » KX*1 

2 KTLLOALMrRI  « *-L  

- SUGGGRIPTG^FOR  PRIMARY  ANGHOR 


V(l»I)  * BBl(l) 


IN  IT  IAL  VALUES 
LEAP  = 1 


PSI  » STAPSI 

P4~  *,rLAI59265 

LOOPE  * 0 

L OOP  A » 

ITEST  * 0 


FFZ(L,2» 
FFXT  L *-3 1 
FFY! L»3) 
FF2T  Lr3» 


0. 

-100. 


i d VfHaiTTaai 


CALL  GFORC 


COMPUTE  CABLE  FORCES  RX»RY»RZ— TENSION  T (M,N)  AND  STRESSED  LENGTH! rtf  Nl 
3 DO  7 NN  » 1*3 
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The  Coiqputer  Program 


SUCCESSIVE 


THIS  PROGRAM  BY  SUB  HASH  C PAHUJA  — - 

THIS  PROGRAM  PERFORMS  THE  MOTION  ANALYSIS  OF  A TRI-MUORED  BUOY 

STRUCTURE-  HAV IMG  AN  AUX  ILIARY  CABLE-43E-IWEEM~IWU~C  ABLE  -LEGS 

SOLUTION  IS  BY  THE  METHOD  OF  IMAGINARY  REACTIONS  AND  SUCCESSIVE 

A P PROX I MAT  I UN  S . . „ . __  ..  . 

THIS  PROGRAM  UTLIZES  THE  CONCEPTS  OF  THE  PROGRAM  DATUBA  WRITTEN 

SKOP  ANO  KAPLAN  OF  THE  NAVAL  RESEARCH  LABORATORY  - - - 

NOMENCLATURE  FOR  THE  MAIN  STRUCTURE 

-44  a.  CABLE  INOEX-  

M = STATION  INOEX 

K = OISCR-ETE  -ELEMENT  INOEX — - — — - 

CABLE  SEGMENT  (M,N) 

BLBAR4M.N)  = UNSTRESSEO  LENGTH  - . - - 

8L(M,N»  = STRESSED  LENGTH 

MCLM+NL-a-  WE-LGHT/  EOQ  T 

XTEN(MvN)  * EXTENS IONAL  RIGIDITY 

MU4M.NI  * DRAG  CHARACRESTIC  - — 

RO'.M,N)  = DRAG  COEFFICENT  RATIO 

KMAX 4 M,  N)  = NO.  OF  04SCRETE  ELEMENTS-  IN  SEGMENT  4M,NI 

KTILDAI M»N I = NO.  OF  DISCRETE  ELEMENTS  IN  FIRST  HALF  SEGMENTIM, 

TIN. NT  a.  TENSION — IN- SEGMENT! H. NT 

( ALPHA! M,N),BETA(M,N),  GAMMA ( M, N ) ) * DIRECTION  COSINES 

DISCRETE  ELEMENT(KtM»N> - ....  

SB  AR(  K,  M,  N ) = UNSTRESSED  LENGTH  FROM  STATION  4M,M  TC  ELEMENTIK 

MEIN.M.NI  = WEIGHT-  -OF-  THE- ELEMENT  - - 

MUEIKyMyN)  * DRAG  CHARATER I ST IC 

STAT  lON-LMyN ) 

4X4M,N),Y(M,N),Z(M,N) ) * COORDINATES 

MMAX-4N1  = NO,  OF  ST  AT  IONS -ON -CABLE!  NT 

(AAKNIy  BBl(N)yCCKN))  = ANCHOR  COORDINATES 

HORTZLTMyNJ  -a  HORTiL  DISPLACEMENT  FROM- THE  GRAV1  TY  PQSI  Tl-UN 
HE IGHT ! My  N ) * VERTICAL  DISPLACEMENT  FROM  THE  GRAVITY  POSITION 

E - -ERROR  FUNCTION 

COMPE  * COMP  AR IS  ION  VALUE  FOR  E 

TIECOM  = COMP  AR  I S-ION  VALUE  FOR  THE  FORCE-BALANCE- CONOLT4  ON 
PSI  = ANGLE  OF  ATTACK  OF  CURRENT 

ST  APS  I = STARTING  VALUE  UF  THE  ANGLE  OF  ATTACK  OF-  CURRENT 
OELPSI  = INCREMENT  IN  THE  VALUE  OF  THE  STAPSI 

- ENOP SI  -=■  LAST-  VALUE-  OF  THE  ANGLE  OF  ATTACK  OF  CURRENT — 

REAL  BLy BLBAR,BLT,MUy  MUE ,MUU,MUUE 

CQMMON/C1/ XF,  X4  2 l,  3 ) y YFy  Y4  2 1,  3)  , ZF  , Z 421  y 3 F - — • 

COMMON/C2/FX! 21, 3I,FY! 2I,3  1,FZ(  21, 31 

CGMMON/C3/W!2W3>yWG!2I,3»y  WEI  10,21,11  

C0MM0N/C4/MMAX ( 3),KMAX( 2 1, 3 ) , K T ILDA ( 21 , 3 ) 

COMMON/C f / BLBAR ! 21y44-y&L4-2-ly3I.  SBARIlQyEl  ,3)  -,1121  ,3!  ,BL  T (2  1 ,31- 
C0MM0N/C6/ AA1(3),BB1(3),CC1( 3 ) , E , DE L T* . JUMP , LOOPE ,LUOPA 

COMMON/C 7/HOR  IZL  (21,3) , HE  IGHT ! 2 ly  3)  - 

C0MM0N/C8/CXC  21, 3),CYI 2l,3),CZ(2I,3) 

COMMON/C9/AV!5»yBV4  5»yVFyV!5),HFyHI5J  ■—  - 

COMMON/C 10/COMPEyCUMPD,  PS  I,  STAPSI ,DELP SI , ENDPSI 

C OHM  ON/ C 1I/-XT  €N4  2Iy-3  > .MU4-2 1,34,  HUE  4 IQy  2Ly  3 Fy  RO  U 1-^34 

COMMON/C  12/ ALPHA! 21, 3), BETA! 21, 3) , GAMMA ( 21 , 3) 

COMMON/ C13/RX(21«3),RY(21,3)yRZ421,i) 

C0MM0N/C14/X0I  21,3),  Y0(  2l,3),ZOI21, 3) 

COMMON/ Cl  5/ XB( 2 U 3 1, YB ( 2 1,  3 ) , ZB  4 21 , 31 

COMMON/ C 16/  FXP(  3 ) , F YP  I 3 ) , F Z P ( 3)  , XP { 3 ) , YP  I 3 1 , ZP  1 3 ) 

COMNON/CIT/DELTAly  POLL  TA  - - 

COMMON/ C 18/ FFX( 21, J ) , FF Y( 2 1 , 3 » ,FF Z ( 21 , 31 


N>  I 

-J 


j 

,M,N»- 


GO  T0(5,4,4)yN 

4 RX4NX,N)  = FXIMX+J44 

RY ( NX , N ) = FY(  MX  , N I 

R7IHX . N>  - FZ4MX.N)  - 

GO  TO  6 

5 RX(NX^i)  ~ RX4MX,  l *+R*V2,24  *RX4  2,  34 
RYlMXtl)  = FY( MX  » 1 )+RY (2*21 +RY 12*3) 

- -RllMXyl)  — — FZ4MX,  i )*KZ4  2,0)  *RZ  40,3) 

6 T4  MX , N)  = SQRT(RX(MX,Nl**2+RY(MX,N)«‘*2+RZ4NXyN>**2) 


8 RX(M,N)  » FX(M,N  )+RX(M+l,N) 

-RY  4N,N)  - a -FY4-M,  N 1 * RYOi*3,N7-  

RZ ( M, N ) a FZ ( M « N ) ♦ RZ(M+ltNI 

GO  TO  16  

220  RX  4 M , N ) = FX  4 M , N ) ♦ RX(M+1,N)  ♦ FFX4L,N) 


RZ ( M * N)  a FZ4M,N)+RZ(M+lyNI+FFZ4LyN) 

14.  = -SQRIIRX  ( My  M)  **2*&Y4M,N  )**  2*RZ  4M,-N)  l 

BL(MyN)  » BLBAR(M,N)*( l.+T(M,N)/XTEN4M,N) ) 


7 BVT<M,N)  = OL4*,N)/T4M,IU 


-DO-40-N  a-  ly4  

MX  = MMAX(N) 

qq  9 H — 2,-t  MX 

X(M,N)  » BLT4M,N>*RX(M,N)+X4M-lyN) 


27  00  110  NN  a 2,3 
-M  1 , NN)  *-X  (NX-r  1 ) 
Y ( ly  NN)  a Y ( NX  y 1) 
1 1 Q-  7 11.  KIM  I 7 ( HX  ■ 1 ) 

X * X Y A f fill  F X V I tN  f 'IT 

10  CONTINUE 


COMPUTE  ERROR  FUNCTION 


LOOPE  = LOOP  E+ 1 


DO  11  N « 2,3 

M s NMAX(N)  — 

11  E * E+(  ( A A 14  N )— X4M ,N I )**2+( BB 14  N ) — Y4  M ,N) I **2  + 4CC l ( N) 
IF(i.GT,COMPE)GO  104  T9,  50, 45),4£AP  — 


- 00  300  N = 1,0  

MX  * MMAX(N) 

00  300  M a 2, MX 

ALPHA! M, N ) = I X( M, N )-X ( M- 1, N ) ) /BL ( M , N ) 

R££A4M,N4  * ( YUM44=-Y-UU--UM  ) ) /&L  ( M,N4  - 

300  GAMMA(M,N)  * ( Z ( M,N  )-Z ( M- 1, N ) I /BL ( M , N ) 


■MM 


ERROR  FUNCTION  COMPARISON  SATISFIED 

— CO 

PRINT  AND  STORE  EQUILIBRIUM  POSITION 


MTEST  * MTEST  + 1 
POELTA  ~~  DELTA 
OELTA  * 0ELTA1 

LEAP  = t 

I F I J TEST  *NE • 1 )G0  TO  3 


JUMP  = 2 

■ - - LOOPS  ■*  0 

00  53  N * 1,3 

MX  = MM  AX  (NT 

DO  53  M * 2, MX 

X04M,N)  -s-  XIM-rNI— 

YOCM.NI  * Y(M,N» 
ZO(M,NI  = ZIMrML 
XB(M,N ) * X(M,N> 
Y«4M»NI  =Y<M*NI 
ZB(M,N)  * Z I M, N ) 

SI  fflMT  tuuc 

GO  TO  61 


COMPARE  ACCURACY  OF  COORDINATES 
52  DO  55  N = 1,3 
DO  55  M = 2, MX 

I-F4-ABS4  X 4MrM>-X04-MrM  1 4-»  GT  »COMPD  »0R  . 

1ABSI YtM,N)-YOIM,N)).GT.COMPD.OR. 

2 ABS4  ZEM,  Nl— ZOIMrN  1 1 .6T  .COMPDIGO — TO-  57 

55  CONTINUE 

ACCURACY  SATISFIED-PRINT  EQQUIL IBR I UM  POSITION 
ITEST  * ITEST ♦ 1 

CALL  T I ELEGIXIL , 2 I , V 4L , 2 I , l < L~»  2 1 , XT L ,31  ,Y(L,3)  ,Z < L ,3  H 
IFILTEST.NE.200IG0  TO  197 


00<-S4--M  - 1-,  3 - 

MX  * MMAX(N) 

00  54  M * 2*MX—  --  

HORIZLI  M,N)  * SQRTI  I XI  M,N)-XB(M,N  H**2  + ( Y (M,N) -YB  (M,NM  **21 

H6IOHTIM,  N>  ^-^4M^N  W-MLH,N1 

56  CONTINUE 

CALL  OVMPOS - - ■ - 

LOOPS  - O - 

LOOP A * 0 


GO  TO  60 

ACCURACY  MOT  AOEQUAT E-RE  ITERATE 

57  DO  59  N » i»3 
MX  » MV  ALU  El  N ) 

— 00  50  -M-»  2+MX  _ 

XOIM, N ) » X I M, N ) 


1 


YOIM*N)  = Y I M*  N ) 

ZQIM.N)  - J1IUW1 

59  CONTINUE 
--••  CO-IO-6I 


* 


C 

c 

c 


ERROR  FUNCTION  COMP AR I SION  NOT  SATISFIEO- 


50  IFI E.LT  «EP IGO  TO  20  - 

INCREASE  IN  ERROR  FUNCTION 
DELTA-—  -DELIA/  2, 


C 

C 

c 


COMPUTE  ERROR-FUNCTION 


12 


13 


OE  * DEL  T A/  SORT!  EP  1 
DO  13  N * 2»3 

MX  - MM  AX  I tl  * 

FXI MX*  N ) = FXPIN)+(AAl(N|-XP(Nl )*DE 
£Y4NX-*N)  = F¥P(NFM~Bai<NI-YPIN14*DE- 
FZ  ( MX*  N ) = FZPIN)+(CCltN)-ZP(Nl  )*DE 
CONTINUE  - 


C 

C 

C 


■H 


CHECK  CHANCES-  IN-  IMAGINARY  REACTIONS 


-DO  14  -N-  =2*  3—  

MX  = MMAXIN) 

— IFIFX4  MX-*N ) .NE  «FXP-4-N  >.QR.- 
1FY(MX,N) .NE.FYPINI .OR. 

2F  l 4 MX  * NT  . NE  . F-Z  PI-N I ICO  TO-  3- 


J 


14 


15 


CONTINUE 

LEAP^  3-  

GO  TO  3 

NO  CHANCE  TIME 
CALL  EX  ITT 

CO -TO-  -LOO- 

DECREASE  IN 

19  LEAP  = 2 

20  EP  = E 
DO  21  N — 2 
MX  = MMAXIN) 


- 

H 


TO  QUIF 


■ 

ERROR  FUNCTUON 

1 

r 3 

A 

21 


X-PTNL- — X I MX  » NT-  - 
YPINI  = YIMX.N) 
ZP4NI  = Z4MX,N) 
FXP(N)  = FXIMX.N) 
FYP4NI  = FY(MX»NI 
FZPINI  = FZIMX.N) 
CONTINUE 


GO  TO  12 


-I 


C 

c 

c 


INCREASE  CURRENT  ANGLE 


60 


61 


19? 


PSI  * PSI+DELPSI 
IFIP&I.r.F  FNDPS  I IGO  TO  » nn 
COSPSI  =*  COSt PSI*PI/180.) 

sinpsi  « siN(Psi*pi/iao.)- 

GO  TO  62 
JUMP  * 2 
LOOPE  * 0 

LOOP  A * 0 

GO  TO  61 


— 4 


62-E 


J 


62  DELTA  « DELTA1 

E AP  — L - 

LOOPA  = LOOPA+1 

6. 

C DRAG  COEFFICIENTS 

DO  30  N * I*  3 

MW  . MM  Kim  

TIT  w I II  > WTTrwr 

DO  30  M = 2.  MX 

4A0GLT«  SORT* ! -BETA!  M,N  TACO SP  SI» ALPHA  I »tN>A4lNP SI  I** 2 

I ♦ GAMMA! M«  N 1**2  I 

BUFFER  * ALPHA!  H,  N) *GOSP^I  ♦ BETA!  MvHIASINPSI 

CX!M,NI=MU!H,NJ*!CAOELT* 

1  L 4-GAMRA!  '♦  BETA!RtNTAA2)*CGSPS4 

2 - ALPHA!M,NJ  * BETA!M,N)  * SINPSI  I 

3 RO!  RtN  ) - A— AL-PHA-TWrN-T- A BUFFER  I 

CYIM,N)*MU(M,NI*(CADELT* 

1 — ( ! GAMMA!  M»N  >»»2  ♦ ALPHALH»H)  **2+^STNPSI — 

2 - ALPHA!  Mt  N I * BETA!M,N>  * C0SPS1  » 

3 - - + R-DLM-t N )-  A -RE TA4H-»N t-  -A  BUFFER  4 

30  CZ!M,NI=MU!M, N l*GAMMAI M, N )* BUFFER* ( RD ! M ,NI-CADELTI 


C COMPUTE  CABLE  DRAG  FORCES  HX» HYf HZ- AND  ELEMENT  DRAG  FORCES  HXE,HYE 
DO  40  N = It  3 

-MX  * RVALUE! NT  - - — - 

00  40  M * 2»  MX 

--  _ — HXHN«MX-N*N 

IFIMXMN-2123,24,24 

24  AT  = AREATMtNt  1 > 

A2  * AREA(M*ltNt  ?) 

NX  -*  GX!H*N)  * A I -A-  GX4**WNI  A--A2 

HY  * CYiM.NI  * A1  ♦ CY ! M* It  N ) * A2 

- HI  « CZIMtNT — *-  AI  &2!~M»Lf-M-T  *-  A3 

DT  * AREA!M»Nt  3 > *ARE  A!  M + l,  N , 4! 

- HX&  *--OT-*  GQSPSI - 

HVE  * DT  * SINPSI 

• GO  TU  26  — 

23  A1  » AREA! MX* It  1 1 

~A2— *-  AREA!2?-2-t-24 — 

A3  * AREA! 2» 3, 2) 

HX  * GXIMtNT  -*  A 1 ♦ CX!  2*  2 I * A 2 ♦ GX4  2t-3l- -A-  A3 

HY  « CY(M>N)  * A1  ♦ CY ! 2t  2 ) * A2  ♦ CY(2f3)  * A3 

HZ  « GZ!  Mt N I * A 1 ♦ CZ!2tZ»  A Ai  ^-WWtW-*-AJ 

DT  « AREA!MX,l,3UAREA!2t2t4H-AREA!2t3t4l 

— - - - WX-E  »— ©T— A CQSPS-I 

HYE  * DT  * SINPSI 

C NEW-  TOTAL  FORCES  - 

26  FX!M,N)  = HX  ♦ HXE 

FYTR,N»  * HY  ♦ HVE  

FZIMVN)  = HZ  ♦ W! Mt  N ) 

40  CONTINUE 

GO  TO  3 

100  STOP-  - . ..  .... 

END 


SUBROUTINC-TLIELEIMXlFyT^ZSFXSrYZ^ZZi 

SUBROUTINE  TIELE6  IS  USED  TQ  SATISFY  THE  FOLLOWING  CONSTRAINTS 

THCCEOMATRICAL  COMRAT  LBTCT XY- 

THE  FORCE  BALANCE 

COORD INAJES  XlrY  WZ  It X 2, Y 2,Z 2 ARE-  TRANSFERRED  FROM -MAIN  -TO -MEET 
THE  FIRST  CONDITION 

THE  SECOND  CONDITION  I S SAT  I SFIEOUS4NO  THE  SEARCH— RQUT-INE 

NOMENCLATURE  FOT  THE  TIELEG  ARRAY 


■ 'Ill- 


CABLE  SEGMENTIM) 

BIN)  * STRESSED  LENGTH 
WWCIMI  = WEIGHT /FOOT  - — 

XXTEN(M)  * EXTENS IONAL  RIGIDITY 


RRDIMI  = DRAG  COEFFICENT  RATIO 

UMAX (X)  NO.  OF  DISCRETE  ELEMENTSIN  SEGMENT !M) 

JTILDA(M)  = NO.  OF  DISCRETE  ELEMENTS  IN  FIRST  HALF  SEGMENT ( M) 

ST4-M-I  = TENSION  IN  SEGMENT  4-MI  

I ALPHA!  M 1 1 BET  A(  M ) » GAMMA!  M ) ) = DIRECTION  COSINES 

-DISCRETE  ELEMENT!  J , M ) 

SAR!  J»M  I = UNSTRESSED  LENGTH  FROM  STATION ! M)  TO  ELEMENTU.M) 

WEE! J, Ml  -WEIGHT  OF  IHE-ELEMENT- 

MUUE 1 J t M ) = DRAG  CHARATERI STIC 

STATION  TNI  — 

(X!M)»Y!M|tZ!MI ) * COORDINATES 

MMAXN  = NO-.  OF  STATIONS-ON  THE  T-IEIEG  ARRAY 

NX  ! 1 If  Y ! L ), Z ( II)  = END  COORDINATES  FOR  THE  TIELEG  ARRAY 
(XtMXIf  Y4MXItZLMX-l-l  — END  COORDINATES  -FOR  - THE  T I E-LEGARRAV- — - 
HORIZLIM)  = HORIZL  DISPLACEMENT  FROM  THE  GRAVITY  POSITION 
HEIGHTIM)  *-  VERT  ICAL  DISPLACEMENT  -FROM  -THE  GRAVITY  -P0S4T1ON- 
REAL  BL • BL  BAR , BLT , MU  t MUE  » MUU,  MUUE 

— DI-MENS  ION  XQ!  224-rYO!22l-rZOI22  > t XRI-32-1  ♦ VS<32 1 , ZB  ! 221  , 

IXP! DtYP! lltZPI 1I»FXP( I),FYP(ll*FZP! 1),AA(22) »BB!22I »CC!22) , 

2W!  22I-»PREMT  22 )»  OISCRTI  221-  

COMMON/H L/FX! 22 1 , FY ( 22  I , FZ 1 22 1 , RX! 22 1 ,RY4 22  I ,RZ (22  I 

COMMON/H2/RX  I!  22  ),RY  I!  22  )»RZ  1(221 

COMMON/H3/ BAR ( 22 1 , B( 22  I , BT( 22 N ST! 221 ,XXTEN(22I 

COMMON/HA/OMAX!  SZT-rST-I  L-OA4-2-2  1 1 SAR  ! 5,221 

COMMON/H 5/ MUU! 221 , MUUE! 5, 22  If  WEE! 5»22l « WWC !22 I , RRD ( 22 ) 

COMMON/H6/ ALPHA! 22lf BETA! 22  I* GAMMAIS2J 

COMMON/ H7/CX! 22I»CY ( 22 ),CZ(  22) 

COMMON/ K8/HOR I ZL( 22), HE IGHT! 221  — 

COMMON/H9/X(22>, Y(22), 2!  22) 

- COMMON/ HI O/F . DEL  TA.I  QOpE .LOOP* 

COMMON/HI l/MMAXN 

COMMON/H12/UUMP, PO I ST,D I ST  - - - 

COMMON/ Cl 8/ FFX(21f3)fFFY(2I,3),FFZ(2If3) 

COMMON/ C 19/ TIECOM.L  

C0MM0N/C9/AV! 5 1 f BV( 5 ) , VF , V(  5 ) • HF , H ( 5) 


rii,riiiTrfAriT/aiiTV|i; 


■ ■ ■ 


I’iM'Ri  I 


COMMON/ C20/ IT EST, JTEST,L TEST, M TEST, L TIE ,MCON 

DATA  9 1/3.141  5924/  

X!  1)  - XI 

VIII  = VI  

Z ( 1 > * Z1 

MX  * -MMAXN  

AA(MX)  = X 2 


non!  . n <■>  o i o r»  n • ! i i f»o 


BBIMX)  = V 2 

cc( mxi  -=  a 

OIST  * SORT ( I AA( MX 1 — X < 1 1 )**2*(BB(MX)-YI 1) 1**2+ (CC I MX ) -Z ( l) I **2) 

HEST  .&€  .1  ICO-fO-^L 

IFIMTEST  .GE.DGO  TO  195 

-BGMPE-*  GOMPEY  2 . 

POIST  * OIST 

C COMPUTE  LENGTH  OF  SEGEMENT 

— _C- \ 

MX  * MMAXN 

TM lyii  jt  My—  i _ ... . . . _ . . __  _ 

00  200  M * 2»MX 

BARI  Ml  » OIST/TMAXN  - j 

200  CONTINUE 


COMPUTE  WEIGHT/ FOOT  FOR  EACH  SEGMENT 


00  116  M = 2?  MX 

PREM4M1  -*  O.  

JX  » J MAX (Ml 

IF<  TQ  ill 

00  115  J = 1,JX 

PREMIM4  * PBEM4-M4-  *-  WEECJ.M1  

115  CONTINUE 

WWC1M1  * PREMI  Ml/ BARI  Ml-  

MWCIM)  = -0.99  * WWC(M) 

nn  r n t|4 

117  WWC(M)  » -0.500 

HE  -CONTINUE-  - 

- COMPUTE  -01-STANCE-  OE  A - O I SCR E TE  SCENE NT--FRQN-  THE  - STA T ION--  _ 


_ 0O-71-N-«-2,~MX 

JX  * JMAX (Ml 

IFTUUEO.OIGO  TO  T1  - 

SUME  * 0.0 

OISCRT  ( M ) * BAR ( M l / ( JX ♦ 1 ) - - 

00  71  K * ltJX 

_ CliMC  f cyiit a.  n | crp  T I M I - 

SARI  Kt  M ) * SUME 

71  CONTINUE  - 

COMPUTE  MIDSEOEMENT  OISCRETC  ELEMENT  J 14  LOATM1  - 


MX  *-  MMAXN- 

00  10  M * 2, MX 
XX  * JMAX4M) 

DO  11  K 3 1»KX 

IF1 S AR4KtMT.GT  .BAR4M1/2.  )GO  TO  10 
11  CONTINUE 

ft  - kx-h __ 

10  JTIL0A1M)  * K— 1 


COMPUTE  STATION  GRAVITY  FORCES  WIM)  AND  INITIAL  FORCES  FX,FY,F2 


MX  * MMAXN 
00  104- M -=  l-.MX 
FXIM1  = 0. 
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106  FVI M I = 0. 

F741>  a 0.  -•  — 

MT  = 0. 

HX  = MMAXN— 1 - - 

00  6 M * 2*  MX 

NX  > 0.  - — - 

KB  * JT IL0A( H ) ♦ 1 
-XX  * JMAX4M) 

IFIKB.GT .KX )GO  10  77 

— 00  - 1^-JC-aL  KB.»KX 

13  MX  = MX+WEECK.M) 

77  JX  » U7lLOA<M*U 

IF4JX.EQ.0IG0  TO  76 
% 00  2 K =l,JX 

2 MX  * MX*WEE( K t M I 
— • — 76 - M4IU-  *~*X  ♦ 4 MMCUX)>BAX4M 
FZ  4 M I a W4  M I 

-6  MT  -=  -MI  ♦ MOO  - -- 

MX  = MMAXN 

FZ 4MXJ  = 1-1*11  * MT 
C 

- C -4NW7AL4ZAI IQU  

C 

SQ2  a SQRT42.4  - 
DELTA1  a ABS4  SQ2*MT I 

- LCON-* .-0--  - - 

LT IE  * 0 

MCOH  -a-  1 

195  DELTA  * OfcLTAl 

JUMP  * l - - - 

LEAP  a 1 

LOOP  A--*-  0 - 

LOOPE  « 0 

■ C COMPUTE  CABLE  FORCES  RX,RY,R Z 


1 

1 

I 

I 

i 

i 

! 

TENSION  TtM,N)  AND  STRESSED  LENGTH 


3 MX  • MMAXN 

RX4MX1  * FX4MX4  - - - 

RYIMX)  « FY(MX) 

PZ1MX1  ~ FZ  I MX  ) 

ST4MX)  « SQRT4RX4MXI**2+RV4MX)**2*RZ4MX)**2) 

B(  MX  I » BAR4MX>*4  1 . *ST 4 MX > / XXTiNFNXl I 

BT4  MX  I a BIMXI/STIMX) 

MX  » MMAXN- 1 - 


00  7 MM  * 1,MX 
- M~»  MX— MM*i  - 


GO  T0I9*  5) » M 
5 RXIM)  « FXIM ) ♦RX4M* 1 ) 

RY4MI  * FY4MI  *■  RY4M  + 11 

RZ  4M J a FZ  4M ) +RZ 4 M* 1 ) 

ST4MI  » SORT  1 RX 4 M ) •♦2+RY4M I **2*RZ  4M| **2 1 
. ami  - nkiiiu  ia<  i ,*cti  m i /kutcmim  i > 

BT4MI  a B4MI/ST4M) 

9 RX4MI  * FX4M|*RX4M*-1 1 

RY4M)  * FY4M1+RY4M+L1 

RZ4MI  a FZ  4 M ) *-RZ 4M*  1 ) 

7 CONTINUE 

COMPUTE  X,Y,Z  COORDINATES  OF  EACH  STATION 
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. - ‘Miiiiiiii  ii  — 


MV  » . MMA^tiil 

I IfV  • II  I ' *■  

IFCE.GT .COHPE I GO  T0(  19, 50, 151 .LEAP 


UP0ATE0  DIRECTION  COSINES 


NX  » MMAXN 

- 00  300  M *-  2tMX  - 

ALPHA!  M I = (X(M|-X(M-1)  )/8(MI 

BE1A4  M4  =-<-Y4-M|"Y(H-lL»/BIM  I- 

GAMMAIN ) = (Z(M)-Z(M-l)  )/B(M) 
300  CONT  INUE 


ERROR  FUNCTION  C0MPAR1  SION  SATISFIED 


51  MX  = MMAXN 

00  335  14 -«  IrHK 
RXl(M)  = RX<  M I 
MY14M1  * RYLNf- 
RZliMI  = RZIM) 


RX1(  1)  * -RX1I1I 
RAF1434  ■«— RV-14  W 
Rzim  = -Rzim 


CALL  SEARCH 

FORCE  BALANCE  NOTORT  AIMED-  START  -RCA  1 N 

IFILTEST.NE.IOOIGO  TO  100 

FORCE  BALANCE- OBTAINEO-PRINT  AND  STORE  EQOI -LIBRIUM  POSITION 


53 


XB4WI  ~X4NI 
YBI Ml  « Y( M I 
ZB4M)  * Z CMI 
CONTINUE 


JTCCT  m ITCCT  a, 

LTIE  « LTIE  ♦ 1 
CAU.  STATS 
MX  » MMAXN 
0053  W * 2, MX 
XOIMI  * XI M I 
V04MI  » YIM) 


CHECK  FOR  FORCE  BALANCE  UNDER  GRAVITY  FORCE S-USE  BINARY  SEARCH  ROUTINE 


LOOPE  = 0 


COMPARE  ACCURACY  OF  COORDINATES 


52  MX  = MMAXN-l 

DO  55  H-  =-  Z*MX — -- 

I F ( ABS ( X ( M I— XO( M ) ) . GT .COMPD .OR . 

1ABS4Y44U-¥04M)4.GI.CGMPO,OR. 

2 ABS  I ZIMI— ZOCMH.GT  .COMPD ) GO  TO  57 
55  CONTINUE 


C ACCURACY  SATISFIED-  PRINT  E QUlLIBRLUN-PO  SI-JION--FOR  -HYDRODYNAMIC-FURC 


DO  38  M = 1*  MX 

RX 1 1 M ) s R x I M I 

RYl(M)  * RY(M) 
RZ1<M>  » R 2-4-M) 

38  CONTINUE 
- RXU  U = -RXU-U 
RYU  U * -RY1(  1) 
R214-I4  = =R2U14 — 


CHECK  FOR  FORCE  BALANCE  UNDER  -AC I-INC 


. GALL  search - - 

FORCE  BALANCE  NOT  OBTAINED- STAR T AGAIN 

FORCE  BALANCE  OBTA  IN-PRINT  EQUILIBRIUM  POSITION  FOR 
MYDERGOYNAMIC  FORCES 
MX  * MMAXN 

DO  5*~M~  - 2*  MX  • - 

H0RI2LI M ) = SORT ( ( X( M I 
HEIGHT  4-M7  — 24H4-2ROO 

56  CONTINUE 
CALL  OYNAMS 
GO  TO  100 

ACCURACY  NOT  ADEQUATE 

57  MX  = MMAXN-l 

00-^7  74  ^-2rHX 

X0( M ) * X I M I 
YG4M1  * VIM) 

ZUIM)  » Z ( M ) 

37  CONTINUE 
GO  TO  62 


RE ITRA7C 


TEST  IF  PRECISION  FOCUS  SHOULD  BE  APPLIEO 


50  IF ( ABS ( AAIMX I-X ( MX  I I .GE.OOMPE.OR.ABS(CCIMX) 

IDOMPEIOO  TO  734  - 

IFCYP(l)  .EO.YIMX  ) )G0  TO  734 

OEL-FYR -4-FY4-HXI»C¥P444-7/ 4 V I MX  l-YP 4-144  *4-604- 

FYP4  1 ) = FYIMX) 

YP41)  * Y4MXI  ......  

FY4MXI  = FYIMX) ♦DEL  FYP 

GO  TO  3 


C ERROR  FUNCT-ION-  COMPART  S I UN  NO-T 


62-K 


JFIE  .LT.EPIGO  TO  20 
INCREASE  4N— ERROR  FUNCTION 
DELTA  * DELTA/ 2. 

COMPUTE  IN AG  I NARY  REE0T40N S 

DE  * DELTA/SQRTt EP  I 

MV  . MM  A VH 

rm  • TnTmliT  ' -----  — — 

FX( MX  I * FXP( 1 I'M A A( MX I-XP!  1 ) )*DE 
F¥-(NXf  * FYPIILM8&I  HXI~*PI  UI*DE 
FZINX)  * F2PC ll+ICCIMX I-ZPt 1))*DE 


• I’TH  “Li 


1 1'  f T*i  I’FT’V- 


IFIFXINX I.NE.FXPI |) .OR . FY< MX  I .NE .F YPt 1) .OR. 

LFZIHXI.ME.FZPL 1IIGO  TO  3 - — 

NO  CHANGE  IN  FORCE  TIME  TO  QUIT 


LEAR  * 3 

GO  TO  3 

45~  IF4MTEST.CE.4 
FFVI Lf  2)  = 2. 
FFVI L*31  * 2. 
GO  TO  100 

tSO  MCON  -2  

.62  CALL  SEARCH 

- CG-TO  LOO 

OECREASE  IN  E 
LO  LEAP  2 
20  EP  = E 

MX  s-  MMAXM  - 

TTM  nnMAfr 

XP(1)  = X I MX ) 


zpin  = z ( mx i 

FXP4  14  - FXtMX  I 
FYP< 1)  = FYIMX  I 
FZP4  14  = FZIHXI- 
GO  TO  12 


r i • * • k # »T •# 1 1 i ■ f a 


30— TO  400 

* FFY (L  » 2 ) 
-FFY4-L.  31- 


ERROR  FUNCTION 


I >am #TT- 


L'tM  rf  "f 


61  cospsi  * cos(psi*pi/iao.i 

SlNPSl  — 4 INI  PSI*P-l  / IRQ.  I 

62  OELTA  * DELTA  1 

LEAP  - 1 

LOOP A * LOOP  A ♦ 1 


MX  « MMAXN 

DO  31  H - 2,  NX  - 

CAPDEL  » SORT  It  BET  A ( M )*COSP  SI— ALPHA IM)*SINPSI I ♦♦2+GAMMA ( Ml **2 I 

BUFFER  » ALPHA!  HL*GOSP  51 * RE4AIM  I*SINPSI 

CXIHI  > MUUf Ml*! CAPDEL*( GAMMA! M ) **2+BETA(M|**2l *CQSPSI  — 


■ T 4*0  |T?'AHIf*>rliT>  IT  <T*ft  #W®  ITT'! 


CYIMI  > MUUI M ) *| CAPDEL*! GAMMA! M I** 2+ ALPHA I MI**2I*SINPSI- 

CZIMI  * MUUI M l*GAMMAI M )*8UFFER* IRRDIMI— CAPDEL) 

31  CONTINUE  . 

NX  » MMAXN-1 

nn  4 1 m — ■>  m v 


A1  * TAREAIM,  1) 

A2  = TAft€A4N+W2>  - 

HX  * CXI N l*A 1+CX IM+1)*A2 

HY  « CY4MI*A1*CY  (M+4  )*A-2 — 

HZ  = CZ I M ) *A 1 +CZ I M+ 1 )*A2 

— OT-a  T ARrEA4Mv-3~4  ♦ TAR£ALM*4. 41 

HXE  * DT  *COSPS I 


o n 


subroutine  input  - - 

REAL  BL,BL8AR,BLT  ,MU» MUE ,MUU, MUUE 

COMMON/ GA/-W1  24-^34  rWC  1 24,  34,  NE44G,-24,A3 

COMMON/ C4/ MM AX ( 3 ) » KHAX ( 21, 3 I ,KTILDA! 21 ,3) 

COMMON/ G5/BL BAR!  21,  3), SC < 2 1,31tSBA«<  IG, 21  ,31 , T121,3)  ,8LT12  4, 34 
CGMM0N/C6/AA1! 3) , SB  1! 3 ) ,CC 1 ! 3 ), E, DELTA , JUMP , LOOPE , LOOPA 

COMMON/ C 9/ AVI  51,BV!  5),VF,V45),HF,H!  51 

COMMON/ C 10 /COM PE, COMP 0, PS  I, STAPSI ,0ELPSI ,ENDPSI 


C0MM0N/H3/BAR! 221,8! 22), BT! 22 ) , ST! 22 ) , XXTEN !22 ) 

COMMON/ H4/JM  AX-4-22),  JTILOA!  221,  SARI  5,221 

' 

COMMON/H 5/ MUU ( 22lfMUUE«  5f 22 ) tWEEl 5f22) ,WWC 122) ,RRDC22I 
rn MHOM/H il /MM AV M 

vunnwYr  n s tl  v 11 V • • — ...  — — - — — 

•—  c 

COMMON/ C19/T  I ECOM.L' 

C COMPARISON  VALUES  AND  CURRENT  ANGLE  REQUIREMENTS 

WRIT  Ei  6, 10 II  - - 

101  FORMAT! • 1* , 15X, • INPUT  OATA  COMMON  TO, BOTH, THE  TIE  LEG  ANji  THE 

1TRI  MOORED  STRUCTURE*/! 

WR  IT  E(  6,  35  ) 

45  FORMATE  SX^'GOMP  AR 1-SIQN  VALUE  S -A NO-  CURRENT- ANGLE  REQUI  Ri-MENJ S*-/ 1 

READ! 5, 1)C0MPE,C0MPD, STAPSI ,DELPSI , ENOP SI , TI ECOM 

1  FGRMAT!F10.2,5FlGrOI - 

WR ITE( 6, 30 )COMPE»COMPD* STAP SI , OELP SI , ENOP SI , T I ECOM 
30  FORMAT!  15X,* COMPARISON  VALUE  FOR  € ^ •-rE15»6-/T-5X,*GOMPARI  SI  ON  VAL 
1UE  FOR  DISPLACEMENT  =',E15.6/15X* * INITIAL  CURRENT  ANGLE  =',F8.3, 


ANCHOR  POSITIONS 

WRITE  4-6, 401  

40  FORMAT! 10X, 'ANCHOR  POSITION  S'/I 

DO  31  N <1,3  - — — - - 

31  REA0(5,2IAAl!N),BBl!NlfCCl!NI 

2  FORMAT  13F10, 01  - 

C NC.  OF  STATIONS  PER  CABLE 

WRITE! 4, 45  I!  AA  1(  N > ,BB  1 ( N 1 ,GG1(N!,N~1, 34 

45  FORMAT! 2X,'AA1» *,F  10.0, 'BB  1 = ' ,F  10.  0 , 'CC  l- ' ,F 1 0. 0/) 

WRITE! 4, 501  

50  FORMAT 4 10X, 1 NO • OF  STATIONS  PER  CABLE'/) 


3  FORMAT!  3151 

WR  ITE  44,  55 1 4-N,  MMAX4  N 1,N*  1 , 3 1 - --  - - 

55  FORMAT42X, 'CABLE  NO.', II, 15) 

REA0!5r344L  

36  FORMAT! 12) 

. ..  up  itch  tn  i 

1 ' " "■  “■  ~~  in  — ■ —m  I.  ii  — ' — — • 

60  FORMAT!/, 5X, 'PROPERTIES  OF  SEGMENTS  ANO  OISCRETE  ELEMENTS' / I 

00  10  N*l,  3 

WR  ITE! 6, 70 )N 

70  FORMAT!  60X,  ' CABLE  NO.',  1-1/ 1- - - .......  

WR  ITE! 6, 65) 

65  FORMAT!  4X,  'M43X, • STRESSED  -L£N..  ' + 3X,  'WEIGHT/E  T-1».3X,.*.EXT»  R I GTI 
1Y' ,2X, • DRAG  CHR ACT  ER ', IX, 'DRAG  COE.  R AT  I O' , 2 X , • K » , 


3 RE NT  ANGLE  =', F8  .3, • DEGREE S '/ 15X 'COMPARI SION  VALUE  FOR  THE  TIE 

44QINT  » *,€46^5/1  

C 

WRITE!4,102!  - 

102  FORMAT4///15X,'  INPUT  DATA  FOR  THE  TRI  MOORED  STRUCTURE*//) 


LEG 
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2 6X»  * S8AR* 1 9X# 'WEIGHT* , 11X, »MUE') 
- _ MMMaMMAX t N 1 
00  10  M=  2»  MMN 


IKMAXCM,N) 

IS  44IRMATC  5FiO_3,F10.2,I2) 

MUCM»N)=1.94*TDR  AG*CABD I A/2  4 
- RO4M*NT=PORAG/T0RAG 
1= M- 1 


iiwn itiii  f tmi j ij» 


1KMAX  C M«  N ) 

24  FORMAT  < 15»  5£  1 5.6»  15) 

IF  < KMAX (M,N).EQ.O)  GO  TO  10 

KKN*KMAX-CM-,NT 

00  5 K-1,KKN 

- - RFAO! 5.  j’QicnAniif.ti 

20  FORMAT ( 4F10 .2 ) 

MU  EC  K,M,  N 1= 1 _94*DRAGCF*XAREA/2.  

5 WRITE! 6, 25)SBAR(K,M,N),WE< K,MfN),MUE CK,M,N) 

25  FORMAT ( 85X*  3E 1 5 .61  - 

10  CONTINUE 


INPUT  DATA  FOR  TIELEG 


WRITEI6, 103) 

103  FORMAT 4 * i-t # 1 5X » ' INPUT-  DATA  FOR  TIE  -LEGAA/T 

WR IT  E( 6 1 37  )L 

37  FORMAT 4 10X, ' -T I E L-EG-UO  INS -CABLS-^AT  STATION  NO. -'-^43/4 
READ! 5,  74JMMAXN 

74  FORMAT!  134  ...  

WRIT  EC  6, 72 )MMAXN 

72  F0RNAT45X,  * NO*  OF  STATIONS -ON  TME- TT£LEG*i-,4  3 )- 
WRITEC6,  104) 

104  FORMAT!  4X.""-a*-'c*T.RIGIfH 
1 AT  10', 3X.'J*,20X,» WEIGHT', 10X, «MUUE ' //) 

. . mv  x yy  Ay m _ _ . ________ 

DO  17  M = 2, MX 

RE  A0C  5,  1 4 ) X XT  E N C M ) , T TDR  AG  , C A BOA  , P P4XUG  ,UMA*CM) 

14  F0RMAT1 4F10.3, 12) 

HUUI  M ) = 1 _Q4*TT  PR  Afi*f  A BOA  / 24  ■ — 

RR0CM)  = PPDR  AG/TTDR  AG 

I * M-l  - - - 

WRIT  EC  6, 16) I«XXTEN(M)»MUUCM) ,RR0(M) , JMAX(M) 

16  FORMAT! 14, 3E 16.7, I 4)  — - - 

IFCJMAX(M) . EQ . 0 ) GO  TO  17 

LKN  * -JMAX4MI-  

00  18  K * 1, LKN 

REAOC  5»29)WEE4K»M), GKAGCF, XXAKEA  

29  FORMAT C 3F10 .2 ) 

MUUECK»M)  * 1.94*GRAGCF*XXAREA/3.  - 

18  WRITE(6.21)WEE(K,M)tMUUE(K,M) 

- 21  F4JRMAT47SX,2E15^63 

17  CONTINUE 

PROVIOE  VELOCITY  PROFILE  IN  SUBROUTINE  VPROFILE 

CALL  VPROFL 

WR1TEC6,75) 

75  FORM  AT ( 40X, • XXXXXXXXXXXXXXXXXXXXXXX • , Tl , *1  •///) 
RETURN 
ENO 


6 

C 

c 


SUBROUT IHE  SEARCH 


c 

--€ 

c 

c 

c 

€ 

c 


tui  q QimftnnT  imp  iMwm  vp ^ thp 
A NEW  VALUE  OF  THE  FORCES  IS 
THE  TWO  PREVIOUS  VALUES 
THE  CONCEPT  IS  SIMILIAR  TO  THE 
FAST  CONVERGENCE 


CONCERTS  OF  8 1 NARV  SEARCH 

FOUNO  THAT  ALWAYS  LIES  INBETWEEN 


DAMPING  FORCE  AND  AS  SUCH  BRINGS  ABOUT 


77 


C 

C 


- CINE  NS  ION  Z^PTjD-rZYR-lZIrZZRCa^ 

COMMON/ C4/MM AX ( 3 I,  KM  AX (2i*3)»KTILDA(2I»3) 

CQMMON/C 18/FFX4  2 1,  3I,FFV4  21,  3I,FFZ  4 21-rST  

C0MM0N/CI9/T IECOM.L 
COMMON/ HI 1/MM AXN 

C0MM0N/H2/RX  1(  22 ) ,RY  U 22  I ,R  Z 1(  22 I 
C OMMON/  C2  O/ITE  SI,  J4GS  T ,T  TE-S-T  ,M  T EST,-LT4C-rMCON— 

GO  T0( 55, 1561 ,MCON 

55  ZXPI II  = FFX4L.21  

ZYPI II  = FFY( L » 2 I 

IIP  I II  = FEZ  1 L, 2 1 

ZXPI  21  = FFX I L , 3 I 

ZYP4  2 1~* — F-FV-4-L-,  411 

ZZPI 21  = FFZIL, 31 

TEST  FOR  FORCE  BALANCE  

MX  = MM AXN 

I F4 1 FFX I L, 2 l+RX 11 I I I »GT .T I£C0M.OR» 

1 4FFY4L,2)+RYll 1II.GT.TIEC0M.0R. 

2 IFF  Z4L,Z4»RZ  1411 1 .61  . T I EC  ON  160- 

76  I FI  4 FFX  I L,  3 I +RX II  MX  I)  .GT  .T  I ECOM  .OR  . 

1 4 FFY4  L,  3 l+RY-14  MX  1 1 »CT  .T  IEGQM  .OR 

2 I FFZIL,  3KRZHMX  I > .GT.T1EC0M1G0  TO 

IR4LTIE.GE.IIG0  TO  50-  

FORCE  BALANCE  UNDER  GRAVITY  FORCES  OB TAI NEO-RE  TURN  AND  PRINT  RESULTS 

-LTiST  * 400  

GO  TO  100 

FORCE  BALANCE  UNDER  ACTING-  FORGES  GBTAINED-RE TURN  AND  PRINT  RESULTS 
90  LTEST  = 200 

GO  TO  100  - - - 

156  MX  - MMAXN 

RY 14  1-1  - = -FRY  4-L-,  2-1  - 

RY14MX)  = FFY I L , 31 
FORCE  BALANCE  NOT  OBT A INED-USE 
AND  START  ALL  OVER  AGAIN 

77  FF  X4  L , 2 1 = -RXlt  1) 

XNO  = ABSIRYli II  I 

YNO  * ABS4ZYR1 1 14 

HNO  = AMIN1I XNO, YNO I 
ZNO  * ABSIXNO-YNOI 
FFY4L, 2 I = ( ZNO/2.  ♦ HNO 1 

FFZ4L,2  I = -RZ II  II  

FFX4  L, 31  = — RX 1 1 MX  I 

X IMG  — - -ABS4RY  4 4MX I 4 

Y1N0  = ABSIZYPI2I) 

HINO  = AM  IN  II  X 1N0,  Y 1N0  I 
Z 1 NO  = ABSI X1N0-YIN0 I 
FFY4L.31  = -IZ1N0/2.  ♦ H1N01 
FFZIL, 31  * -RZ1IMXI 

100  RETURN  - - ... 

END 
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BINARY  SEARCH  TO-OBT44-N  NEW  VALUES 


SUBROUTINE  GFGRC  

REAL  BL»  BLBAR, BL T, MU, MUE .MUU.MUUE 

COMMON/ G2/ FX I 24r-34rF»4^1f  34tFZI  21,34 

COMMON/ C3/M( 21,3), MCI 2 ly 3 I , WE* 10. 21 , 31 

COMMON  AG4/MNAX 4 3 I.KMAXI  2 U 34»KT  IL0A4E1tBT  — 

COMMON/ C 5/ BL  BAR ( 21. 3) ,BL( 21. 3), SBARI 10,21 .3) .TI21 .3) .BL1 
COMMON/ C4/AAH  3)fBB44  34  ,CC  14  3),  E . OEL  TA  . JUMP -ylQQFE  ,IQGPA 
COMMON/ C 17/0ELTA1.PDELTA 


MT  » 0. 

DO-  3 — N •■■■■If — 3-  ■ — •—•■■■■  — ■ - — - - — — 

MX  = MMAX(N)  - 1 

- 00  3 M-»  2,— MX  — 

UX  * 0. 

f*EMr  N4—  «-<>, 

FYCM.N1  = 0. 

KB  -*  KT  I4.-OA4M.NT  *-  1 --  — 

KX  = KMAXIM.N ) 

1F4KB.BT  *KXI  GO  W W-  

00  1 K * KB.  KX 

10  KX  = KT IL0A4 M* 1. N I 

- IF4*X»§Q.01  GO  TO  30  - - 

00  2 K - 1.  KX 

2 MX  * -MX  ♦ WE4K-.N*  1-,  N I 

30  M(MyN)  - WX+MC(M,NI*BlBAR(M,N)/2.+WCIM+l,NI*BLBAR(M+i,N)/2 
—•  «4MfN4  * M44M4T 

3 MT  » MT  ♦ W( M.N) 

MX  -*  0.  • - - - • 

FXI MMAX4 1), 1)  = 0. 

FY 4 MMAX 4- 1 ) . 14  * -0.  

KB  * KT ILDA{ MMAX ( 1). 1)  ♦ 1 


IF(Kfi.GT.KX)  GO  TO  20 

DO  4-- 44-* • KB,  KX  - — - 

4 MX  = MX  ♦ WE(K,MMAX(  1). 1) 

20  CONTINUE  - 

00  6 N = 2,  3 

KX  - KT  4-L  0A4  -2fN  ♦ 

FX( MMAX ( N I , N I = 0. 

FY  I MMAX4-N-)  . N 1 = 0.  

IF(KX.EQ.O)  GO  TO  6 
00  5 K = 1,  KX 

5 MX  = MX  ♦ ME  ( K,  2. N ) 

A MX  = MX»MC  ( 2 ,N4  *BLBAR42-,N4/-3-« 

M(MMAX(1),1)  = MX*MC(MMAX(l)yl)*BLBAR(MMAX(l) 
FZ (MMAX4  1).  1)  * MI  MMAXI  1>*1T 
DO  35  N * 2,3 

FX4l.NI  = FX4MMAX4 1>. 1»  - 

FYIl.NI  * FYI  MMAX(  11,11 

35-FZ44,N4  - F^IMMAXTHf  IF  

WT  = WT  ♦ WtMMAXI 1),  II 
00  7 N * 2,  3 
7 FZ (MMAX ( N ) • N I = -MT  / 3. 


C 

C 


COMPUTE  INITIAL  DELTA 
SQ2-  « SORT  4 2.  I - 
OELTAl  * ABSI SJ2*WT/3.  I 


DELTA  * DELTA1 

RETURN 

END 
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MICROCOPY  RESOLUTION  TEST  CHART 


t 


CUMClUUll  iflFil  M.  M.  Ififl  > 

• VUV  V IUIV  WfTV  R | M | f ^ f I VU  V —————  — — — — . 

REAL  BL,  BLBAR,  BLT,MU»MUE,MUU,MUUE 

C0MM0N/Cl/*F-rX4ai-r3I*-¥F-,V-l21,-3>tZFtZ!21.3) 

COMMON/ C 4/ HH AX I 3), KM AX! 21. 3 ) .KTIL0A1 21 ,31 

COMMON/ 05/ SL  OAR4-2 1 ,3»,OL!2  lr3  )*-  SBA4U  10^24^3 1 421-^  - 

C0MM0N/C9/  AV(  51 , BV  ( 5 1 , VF,  V<  5 >,  HF  ,H  ( 5) 

4U3MMON/C11/XI  EN<  2 W i * , MU!  2 W 3 WMUE ! 10,21  *31  *RD!21,31 

COMMON/C  12/ ALPHA! 2 1,31, BETA (21, 3 1 .GAMMA ! 21 , 31 

CO- TO!  100,200,300,400)  t4W 

C LINE  INTEGRAL  BELOW  STATION 

- 100  CONTINUE  — 

KL  * LIMITIZI M- 1,N ) +GAMMAI M »N)*BL! M, N) /2« I 

...  . . KU-  * LIMIT  -I — Z4-M»N4  -4 — ■ — . - 

TOP  = BLI M, N ) 

-BOT-*  BUM.NI/2. 

2 KMIN  * MINO  ( KU,  KL  ) 

KF4JP  = -KNIN  -»•  l * I ABS  ! -KlMCO-l 

SUM  * 0. 

1F4XNIN.GT  .KT0PJ--4ia-4O-  lO-  

00  1 K - KMIN,  KTOP 

-XT-*  ! HI  K )— Z I H-  l.N  1 1 / GAMHAIM.N) 

1 SUM  = SUM  ♦ AREA SIM, N,K,XI)  - AREAS! M,N  ,K+1 , XI 1 

RETURN 

C - -LINE-  4NTSCNAL-  ABOVE  -STA-T-ION 

200  CONTINUE 


KU  * LIMITIZtM- l.N ) ♦GAMMA! M ,N)*BL(M,N)/2.) 

4 OP  » *04NrN4A2.. 

BOT  = 0. 

. . go  jq  2 - 

C ORAG  OF  ELEMENTS  BELOW  STATION 

30O--GQNT4NUS 

KB  * KTILOAIM.N)  ♦ 1 

-KX-  * -K  MAX-I M,  N ♦ , 

5 SUM  - 0. 

IFIKB.&T  .KX I GO  TO  40  - 

DO  4 K * KB,  KX 

• - ZT  -»  ZIM-UNI  » GAMMA1M.NI  ■*  SOAR  I K , M,NI«  1 1 .♦  TI M,  N) /XT6NIM.NI  1 

J « LIMIT!  ZT  ) 

4--SUN  »-  SUM— *•-  MUE!K-,M,NI  * 14V4  J4  4 4T4**2 

40  AREA  « SUM 

- RETURN  - - 

C ORAG  OF  ELEMENTS  ABOVE  STATION 

400-  CONTINUE 

KB  « 1 

-KX  * KT  JLDA1N.N ) - 

GO  TO  5 

END — ■-  ...  - 
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FVHCHQM  A*EAS4*MML*X-I  + ~ 

COMMON/C  1/XF.  XI 21.3l.YF.YC 21. 3 1 . ZF , Z C 2 1 .3* 


COMMON/C 12/ ALPHA! 21.3)»BETAI 2i.3l.GAHMAI2l.3l 

niB^^0Ik7H7Nr.  «.  , 

♦LBVilU*GA*MA4*r*i)**^  * 4-Xi««3/3^l 

RETURN 


SUBROUTINE  VPROFL  — 

COMMON/ C9/  AV ( 5 1 » BV(  5 » , VF  , V!  5 ) ,HF  ,H  ( 5) 

1 READ(5'2)H(K),V(K) 

2 FORM  AT! 2F10 .0 ) — - 

IF  ( H(K)  .GE.  40000.  I GO  TO  3 

K * K -♦  I — - • — - — -j 

GO  TO  I 

^ KX~» - K 

DO  4 K=2,KX 

BV4K!  = 4 V4  K )-*V4  K”  II)  / 4HilO-H4X^444 - - 

AVIKI  = V(K-1»-BVIK>*H4K-1I 

4 WRlTE!**9TAV(K>t»VFX>  ' * * • 

9 FORMAT!  10X,6HAV(K  )=  , E16 . 5/ 1 OX, 6HB V 4KI  = ,EI6.5) 

~~  ....  RSTUR* j 

END 


FUNCTION  TAREA <M,HGO> - - 

REAL  BL  » BL  BAR  , BL  T,MU,  MUE  ,MUU, MUUE 

i- COMMON/H3/BARI  22-F,B4224tOTI-224t-ST-4-22)  , XXTEN  42-24 

COMMON/ H4/JMAX 4 22) , JTILDA! 221 . SAKC 5,221 

!-  €OMMQN/H5/MUU4-224,MUUEIS,22T,WEE4  5T224^WWG422),RRt)l22)  - 

C0MM0N/H6/ ALPHA! 22 ) » BETA! 22  ) , GAMMA! 22 I 

COMMON/ H9/X4  224*Y4224'rZ-422  ) - 

COMMON/ C9/ AV!5),BV!5)«VF,V!  5),HF,H!5) 

An  T nl  1 AAA  OAAA.  3Anfi  _ LAAft  I , MAH 

W TT/t  »WVT  C WVf  9 VVVt  “ VW  T f n\|LT’ 

C LINE  INTEGRAL  ABOVE  STATION 

l/jnn  rnki  iuuc  ....  

KL  * LIMIT!  ZIM-|UGAMMA(ML*BJH)/2.J  

KU  « LIHIT4Z4MM  - - 

TOP  * B! M > 

..  flfiT  t -- ■_  . ...  

2 KMIN  * MINOI KU, KL I 

KTOP  * K4IIN— 1*1  A&S4  KU-KL4  

SUM  * 0. 

IF4KMIN.GT.KTOPJGO  TO  10  - — 

00  1 K = KMIN, KTOP 

-XI  - LH4  K )-Z4-N*-4  M/ GAMMA!  M I 

I SUM  = SUM+T AREAS ! M,K , X I I— TAREASIM  »K+ 1 , X I ) 

10  TAREA  = SY6N4KU-KL  l*SUM  ♦ T AREA S!**KU,  TGP4-TARE  AS 4M^KL  ,BGTT 
RETURN 

C LINE  INTEGRAL  BELOW  STATION  — - 

2000  CONTINUE 

....  ...  Mi.  . j£-  l iu  L T i 7 l y — | II  . . , _ . . — . 

KU  = LIMITIZ1M-1I +GAMMAI M|*B!M)/2. ) 

TOP  = B4MI/2.  - 

BOT  = 0. 

GO  TO  2 - - - 

C DRAG  OF  ELEMENTS  ABOVE  STATIONS 

!-  -MOO  CONT4NUE  - — 

KB  * JT  ILOA! M I ♦ l 

KX  * UMAX!  M ) 

5 SUM  = 0. 

TF4KB.GT.KXIGO  TO  20  - 

DO  4 K = K B»KX 

J * LIMIT4ZT) 

4 SUM  * SUM+MUUE4  K,M  )*!AV(4)*8V!  J4*Z  T)**2  

20  TAREA  = SUM 

RETURN  - 

C ORAG  OF  ELEMENTS  BELOW  STATION 

4000  CONT4NUE-  - -- 

KB  * I 

KX  * JT  ILOA 4 M ) - ---  - — 

GO  TO  5 

ENO  - - 


RIHfTHlW  I IMfTi  tn 

C0MM0M/C9/AV( 5>»BV( 5)»VF»V( 5)»HF»H( 51 

flQ-t-  4 - aw§ 

IFIZT— H( J I )2t  2? 1 

- — l ^GWWNUe 

2 LIMIT  « J 

AfilUAAl 

nfcrwnlf  — — “ — — • — — — — ■— 

END 

I 


l7iTrfWlITfIT 


COMMON/ C2/FX! 2 1*  3 ) • FYC 2 1, 3 ) . FZ ( 21 , 31 

COMMON/C 5/ BL BARI 21*  3 ) * BL C 21 , 3 ) * SBAR 1 10, 21 ,3) ,T(2l ,3)  *BLT(21 ,3) 


COMMON/ C7/H0R IZL I 2 1, 3 1 * HE IGHT< 2 1, 31 


COMMON/ C 1 1/ XT EN (21, 31* MU (21*3), MUE ( 10, 21 ,3) ,RD(2 L ,3 > 


COMMON/ C 17/ DEL T A 1,P DELTA 

• * • • « « > • 


5 FORMAT!  IX, 41HEQUIL IBR IUM  POSITION  UNDER  GRAVITY  FORCES ,// 10X , 

2 ,15) 


MX*MMAX!N) 

DO  10  Ma-2«MX 


10X,8HFX!M,N)*,E16.9, 10X , 8HFY I M, N) *,E 16. 9 ,1  OX ,8HF Z I M , N> - ,E 16 .9/ 


639  FORMAT  I 10X, 8HRXI M,N )», E 16.9, 10X , 8HR Y !M,N ) »,E 16.9 , 10X, 


GO  TOC  10, 9), JUMP 

9 WR-TTiCA,  LI  *HORTZC4JLrN4,WCICHTCM^N-l 

11  FORMAT ( 15X, 12HH0R IZL CM,N )*, E 16. 9, 10X, 12HHE IGHTI M,NI  = ,E 16. 9 1 

10  CONTINUE  

RETURN 


WRITE16, 13) 

_ EC  6 , 1 16 ) - 

116  FORMAT!  13X  , • EQU IL IBR  IUM  POSITION  Ul 
17  WR  IT  €4  4 2 »E,  OEL-T  A,LOOP  £^4.  OOPA,RG4 

12  FORMAT! 10X,2HE«, E16.9, 10X,6HDELTA= 


2F8.3/ ) 

GO  TOA 
ENTRY  EX  ITT 
WRITEi*,13) 
WRIT  EC 6, 14) 


1 TO  CHANGE  THE  IMAGINARY  REACTIONS 

2/ LX  *TOO  NEITHER  ACCURACY  REOUIREMEMT-S-ARE-  TOO-  SMALL  C COMP £7-01 
3LE  HAS  GONE  SLACK  CCHECK  TENSIONS). 

4/ IX, 55HPR INTOUT  IS  G7VEN  FOR- TROUBLE  SHOOTING -PURPOSES  ONLY 
GO  TO! 16, 17), JUMP 

1A  FORMAT  I 1H1-4 

END 


1 O O I J t o o o 


REAL  bl,blbar,bl;,mu,mue,muu,muue 

COMMON/ Hl/FX!  22),FV!  22)tFZ4  ag-MW  221  tRY422)  »RZ- 

COMMON/H3/ BAR( 22 ) , B< 22 ) v BT( 22 ) t ST( 22) * XXTEN (22 ) 

COMMON/ H8/  HOR  I ZL  ( 22 ) , HE  IGHTI  22  ) 

COMMON/ HIO/ E* DELTA*  LOOP E ,LOOPA 

COMMON/ Nil/ MM AX N 

COMMON/HI 2/ JUMP »PDIST»DIST 


OUTPUT  -FOR- 


WRIT  E(6» 13) 

13  FORMAT-!- 4W1) 

WRITE! 6t 100 )POI ST 

100  FORM AT4-40Ry,4EN0T-H  -OF— TM6-  —41  EL€G* 1 ,616.9/4 

26  WRITE!6,25)Et  DELTA* LOUPE 

— 25-  FORMAT4-4XA1M6QU1L  IBR  IUM  POSITION  UN08R  GRAVITY  FORCE  S ,// IPX  , 
12HE=,F16.9, 10X, 6H0EL TA=, E 16. 9, 10X, 1 9HN0.  OF  ERROR  LOOPS*«I5l 
... H-  WR1 TE46,  104  M41ST  

101  FORMAT! 1 OX, 'CHORDAL  OISTANCE  BETWEEN  ENO  POINTS*' ,El6.9/i 

. MW  *.  |J  _ _ . _* . ----- 

DO  20  M = 2»MX 

M— L 

WRITEI6, 28) I* FX( M )« FY(M ) ,FZ ( M) , X(M) ,Y(M) ,Z(M) , ST (M) ,B(M) 

28  FORM  AT  4 / » SX-t  ^SFGMEMT  -NUMBER  »«t12^/4Q*-, 

1 6HFX!M)=,E16.9,  10X»6HFY(M)*» E 16.9, 10X*6HFZ(  Ml  «,E  16.9  ,/10X, 

3 6HST(M)*,E16.9, 10X, 5HBI  M) *, E 16. 9) 


1000  FORMAT! 10X* 6HRX ( M )*, E16 .9, 10X, 6HR Y( M ) =, E 16. 9, 10X, 6HR Z (M) * , E16 .9 » 

GO  -TO  T 2Q*-2B4r JUMP- - 

29  WRITE! 6, 11) HOR IZL(M), HEIGHT (M) 

11  FORMAT  ( 15X 1 10HM0R  IZL  4M  > = t E 18.P*  10X-»  1 OHHE-46H-T4M)  » *€  16.  94 

20  CONTINUE 


OUTPUT  FOR  FOU4L4B8-4UM  PO  SI  T ION  -UNOfeR 


ENTRY  OYNAMS - - - - - 

WR  IT  E(  6 . 13) 

_ yft  ITE(6|  16) 

16  FORMAT!  10X,' EQUILIBRIUM  POSITION  UNDER  ACTING  FORCES'//) 

17  WR  IT £46*42)  E*OE4TA*LOOP€*iOOPArF41- 

12  FORMAT! 10X,2HE*,E16. 9, 10X, 6HDEL TA  = » E 16. 9/13X, 1 9HN0.  OF  E RRCR  LOOPS 

l*.  I^»  1QX.22HNO.  OF  ACCURACY  LOOPS=~,l  5Y-LSX*!  4HC  UR  RE  NT- AMCLF^* 

2F8.3/) 

-GO-TQ-37- 

ENO 
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D.  SAMPLE  OUTPUT 


A computer  OUTPUT  as  obtained  from  the  subroutine  OUTPUT  and  TIEOUT  is 
reproduced  in  this  section.  This  is  done  to  present  the  format  in  which  the 
output  should  be  expected. 

>c«  # , « • • • • • • • • .-•••#«  • 
t First  part  consists  of  the  equilibrium  configuration  of  the  array  system 

when  only  gravity  forces  are  acting.  Next  follows  the  configuration  under 

the  action  of  hydrodynamic  forces  produced  because  of  the  current  (in  this 

case  only  zero  degree  angle  of  attack  is  considered).  The  printout  for  main 

cable  array  is  shown  only  for  cable  number  1.  There  is  a similar  format  for 


SEMI*  NT  KUNRERMO 

_ . FM"»-  0.0  ._.  __  rriR|«  0.0  ...  FZI0I-  Q.S683302T3E  01 

XINI-  0.883T134T2E  04  Tl  HI  —0.  48847698  IE-01  ZIM-  0.412474121F  03 

STCNl.  0.1 044338 1.3E  04  8t*l>  0.153133669E  04 

RXIH»«  0.0  RFIRI-  0.104433444E  04  R2INI>  0.284390230E  01 
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APPENDIX  II 


SPECIAL  FEATURES  FOR  USING  THE  COMPUTER  PROGRAM 

The  Computer  Program  that  has  been  described  in  Appendix  I can  be  used 
by  using  appropriate  data  cards.  Because  of  the  truncation  errors  made  by 
the  computer  it  is  necessary  to  choose  the  appropriate  cutoff  values  that 
vill  suggest  that  the  method  has  been  used  satisfactorily  to  simulate  the 
cable  system.  A brief  description  of  these  cutoff  values  is  made  in  the 
following  paragraphs: 

A.  Cutoff  Value  to  Define  the  Acceptable  Completion  of  Imaginary  Reaction 
Routine 

As  described  earlier  in  Chapter  I,  the  equilibrium  configuration  of  the 
array  system  is  obtained  if  E,  the  measure  of  error  between  the  end  coordinates, 
as  obtained  by  system  of  forces  at  hand  and  the  correct  ones  as  specified,  is 
zero.  Theoretically,  the  iteration  would  continue  until  E is  identically  zero 
i.e.,  until  the  equilibrium  configuration  (under  the  constant  applied  forces) 
is  obtained  exactly.  However,  because  of  computer  roundoff  errors,  it  is  not 
possible  to  achieve  this  and  as  such  a cutoff  value,  that  would  define  the 
, acceptable  completion  of  the  imaginary  reaction  routine,  has  to  be  defined. 

If  this  cutoff  value  is  denoted  by  COMPE,  then  the  imaginary  reaction 
routine  is  considered  to  have  given  the  satisfactory  equilibrium  configuration 
of  the  cable  system  when 

E - COMPE  (1,11) 


That  is,  when 


^(nj.n:  YM(n),n:  ZM(n),n  for  n - 2 and  3 for  the 
main  arrays  and  # for  the  tie  leg  array,  are  all  within  COMPE 

from  their  true  anchor  values. 
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B.  Cutoff  Value  to  Define  the  Acceptable  Completion  of  the  Successive 
Approximation  Iterations 

As  described  in  Chapter  I COMPD  denotes  the  fixed  accuracy  value  and  is 
used  to  suggest  the  acceptable  completion  of  the  successive  approximation  routine. 
That  is , the  equilibrium  coordinates  of  any  cable  station  for  two  successive 
iterations  are  compared.  If  the  coordinates  differ  by  less  than  COMPD,  the 
iteration  is  considered  satisfied. 

However,  introduction  of  COMPE  as  a cutoff  value  also  introduces  an  error, 

again  within  /COMPE  from  the  actual  equilibrium  coordinates , into  the  coordinates 

calculated  for  every  cable  station.  With  the  result,  it  is  important  that  the 

cutoff  value  COMPD  for  the  successive  approximation  iteration  be  chosen  outside 

the  value  of  this  inherent  error. 

(9) 

A safe  minimum  value  for  COMPD  is 
COMPD  * 10  /c6M£E 

C.  Cutoff  Value  to  Define  the  Acceptable  Completion  of  the  Force  Balance 


Condition 

The  force  balance  condition,  as  described  in  Chapter  IV,  is  said  to  have 
been  satisfied  if 
FFX  + RX1  = 0 
FFY  + RY1  = 0 
FFZ  + RZ1  * 0 

where  FFY  and  RY1  are  defined  by  equations  (1,4).  FFX,  FFZ  and  RX1,  RZ1  are 
similar  forces  in  the  X and  Z direction  respectively. 

However,  because  of  computer  roundoff  error  some  cutoff  value  is  defined 
which  would  suggest  that  the  force  balance  condition  has  been  satisfied  to  an 
acceptable  degree.  This  value  is  defined  by  TIECOM. 
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TIECOM  can  be  made  very  small,  since  binary  search  routine  will  make  the  values 
of  FFY  and  RY1  converge  to  an  acceptable  positive  number.  An  acceptable  value 


for  TIECOM  could  be  1 lb. 


D.  PRECISION  FOCUS 


This  convergence  concept  was  suggested  by  G.  H.  Savage  'and  is  used 
in  conjunction  with  the  6 method  of  convergence  discussed  in  Chapter  I.  It 
utilizes  the  past  performance  of  the  force  and  displacement  to  bring  the  cable 
array  to  the  required  coordinates.  Referring  to  Figure  17,  the  end  of  the 
cable  is  released  from  point  K'  and  imaginary  reactions  are  applied  at  the 
free  end.  For  each  iteration,  depending  upon  the  position  of  the  K end,  additive 
forces  given  by  equation  (20,l)  are  added  to  the  imaginary  reactions.  After 
a certain  number  of  iterations,  it  is  possible  to  bring  the  K end  to  its  X 
and  Z coordinates  within  limits  of  COMPE.  Normally  the  above  process  would 
be  applied  until  y coordinate  is  also  within  the  limits  of  COMPE.  However,  it 
has  been  found  that  this  situation:  is  not  always  possible  to  achieve.  The 
positive  number  6 often  becomes  too  small  to  make  any  significant  difference 
in  the  additive  force  given  by  equation  (20,l),  the  E never  becomes  less 
than  COMPE  and,  therefore,  equilibrium  configuration  is  not  obtained.  To 
cope  with  this  problem,  after  x and  z coordinates  have  been  obtained  within 
the  limits  of  COMPE,  the  Precision  Focus  convergence  approach  is  applied.  To 
demonstrate  this  method.  Figure  17  represents  the  plot  between  force  FY, 
which  is  the  end  reaction,  and  the  Y coordinate.  Pt.  H represents  a point 
of  intersection,  the  point  which  represents  the  exact  y coordinate  to  which  the 
k end  of  the  cable  should  strive  to  reach. 
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FORCE  FY 


Figure  17 : The  Precision  Focus  Concept 

If,  at  the  present  iteration,  the  y coordinate  (Yp)  to  which  K end  of 
the  cable  has  reached  is  represented  by  point  P on  the  curve,  then  at  this 


point,  the  force  Fy  pulling  at  the  end  is  given  by  OQ.  Also  at  this  point, 
a record  is  kept  of  previous  iteration,  wherein  the  y coordinate  of  the  K 
end  of  the  cable  was  represented  by  point  M and  a force  Fy  represented  by  ON 
was  pulling  on  it. 

Then  between  these  two  iterations,  MP  represents  the  slope  of  the  curve. 
If  6 is  the  angle  as  shown  then 


tan  6 = 1 


where  GP  = OQ  - ON 
and  MG  = MN  - PQ 


Also  at  this  point,  the  aim  is  to  increase  the  y coordinate  by  (yp  + PQ) 
and  to  do  this,  a force  increase  represented  by  HQ  will  be  required  to  do  so. 
Thus,  the  new  additive  force  is  given  by 

- § * k> 


This  force  is  added  to  force  (Fy)  which  coupled  with  forces  FX  and  FZ 
would  bring  the  cable  to  point  R and  the  process  is  repeated  until  y coordinate 
is  obtained  within  the  limits  of  COMPE. 
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Since  the  above  approach  utilizes  the  past  performance  and  takes  into 
account  the  slope  of  the  curve,  this  method  seems  to  avoid  the  problem  posed 


by  the  positive  number  6 getting  too  small. 

Close  to  point  H,  the  curve  gets  flat  with  the  result  that  it  is 
possible  to  get  into  a position  wherein  points  P and  M are  same  i.e.,  PQ  - MN  = 0,  * 
When  this  happens  the  control  is  reverted  back  to  th-  irevious  convergence  concept.; 
This  concept  is  utilized  until  x and  z coordinates  are  within  the  limits  of 
COMPE,  at  which  point  control  is  passed  over  to  Precision  Focus  again.  This 
provess  is  then  repeated  until  E <_  COMPE  at  which  point  equilibrium  configuration 
is  obtained. 

E.  ASSUMPTIONS  USED  IN  THE  COMPUTER  MODEL* 

In  this  section  a discussion  is  made  of  the  assumptions  made  in  the 
computer  model.  These  assumptions  are  made  only  for  the  purpose  of  convenience 
and  are  not  necessary  from  a theoretical  point  of  view. 

1.  Subroutine  VPROFL  deals  with  the  current  profile.  In  the  program, 
a profile  assuming  a series  of  straight  line  segments  of  arbitrary 
length  and  scope  has  been  used.  However,  the  analysis  developed  has 
been  unrestricted.  The  analysis  requires  evaluation  of  the  integral 
as  defined  by  equations  (29,3).  Thus,  any  velocity  profile  could  be 
used  as  long  as  the  above  integral  could  be  evaluated. 

However,  if  current  profile  consists  of  series  of  straight  line 
segments,  as  in  our  case,  this  integral  can  be  evaluated  exactly**. 

As  such,  this  profile  was  chosen  to  write  the  computer  program. 

2.  On  all  runs  made  with  the  computer  program  as  given  in  Appendix  I,  a 
uniform  current  profile  was  used.  Current  velocity  of  0.6739  per  sec. 
was  assumed  all  the  way  from  the  bottom  (z  = o)  to  the  surface.  This 

*The  reader  who  desires  to  use  the  computer  program  using  a different  data 
set  than  listed  in  Appendix  I is  referred  to  Sub-Sections  6 and  7 of  this  Sectifl 
**The  evaluation  is  performed  in  Appendix  II,  Skop  and  Kaplan,  "Static  Configuratl 
of  a Tri-Moored,  Subsurface  Buoy." 
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velocity  profile,  as  suggested  earlier,  was  broken  up  into  series  of  straight 
line  segments.  The  program  is  written  to  take  into  consideration  any  velocity 
profile,  but  a constant  profile  was  used  in  this  study  to  provide  easy 
comparison  of  hydrodynamic  forces  at  all  points  on  the  array  system. 

3.  Only  one  angle  of  attack  of  current  was  used  to  study  the  hydrodynamic 

behavior  of  the  structure.  The  facility  to  have  & different  angle  of  attack 
exists  in  the  program. 

It  was  observed  from  various  computer  runs  made  with  the  computer  program 
as  listed  in  Reference  9,  and  modified  to  suit  the  facility  at  the 
University  of  New  Hampshire  .that  -the  wrst  displacements  of  the  Subsurface 
Buoy  were  encountered  with  the  angle  of  attack  of  the  current  as  zero  degrees 
(6=0).  As  such,  for  the  analysis  of  the  present  array  system,  current  with 
only  this  angle  of  attack  was  used. 

U.  As  described  in  Chapter  IV,  a Binary  Search  Routine  was  suggested  to  be  used 
to  satisfy  the  force  balance  condition  in  aminimum  number  of  subroutine 
iterations.  Forces  in  the  y - coordinate  directions  posed  the  majority  of 
problems.  As  such,  this  routine  was  used  only  in  the  y - direction.  Should 
continued  use  of  this  program  indicate  difficulties  in  the  X or  Z coordinate 
directions,  then  the  Binary  Search  Routine  may  need  to  be  applied  there  too. 

5.  The  results  that  are  discussed  in  Chapter  V were  obtained  by  using  some  of 

the  data  listed  in  a report  by  R.  G.  Paquette  entitled  "Seaspider  Hydrodynamics," 
dated  2k  April  1969.  This  uses  a surface  buoy  along  with  the  subsurface  buoy 
in  accordance  with  the  lumped  parameter  representation.  This,  however,  is 
not  necessary  and  any  other  data  set  can  be  used. 

6.  The  length  of  the  tie  leg  array  and  length  of  each  segment  of  the  tie  leg 
array  were  determined  by  the  program  internally.  This  also  was  done  only  for 
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simple  convenience.  If  the  length  of  the  tie  leg  is  to  be  given, 
then  a change  in  the  program  would  have  to  be  made  accordingly. 

7.  The  tie  leg  array  is  used  as  a neutrally  buoyant  array.  However, 
since  it  is  not  possible  to  have  this  in  practice,  it  was  decided 
to  make  the  array  1%  positively  buoyant.  This  was  done  by  using 
floating  objects,  with  their  weights  in  water  as  a reference. 

From  this  the  weight  of  cable/ft.  in  water  was  determined  internally 
by  the  program.  This  situation  can  very  easily  be  reversed.  Also, 
since  the  length  of  each  segment  was  determined  internally,  it  was 
also  decided  to  space  the  floating  objects  equally  in  the  segments 
by  the  program  internally.  This  is  completely  for  the  purpose 
of  convenience, 

F.  OTHER  CONSIDERATIONS 

Besides  the  above  cutoff  values  there  are  special  cases  in  which  the 
program  may  never  obtain  the  equilibrium  configuration  i.e.,  E may  never 
obtain  the  value  of  COMPE.  There  can  be  one  of  the  two  possible  reasons 
for  this  to  happen. 

1.  The  convergence  factor  6 , as  discussed  in  Chapter  1,  approaches 

zero  as  E approaches  zero.  Consequently,  for  a very  small  value  of  6 
no  change  will  occur  in  the  applied  imaginary  reactions  as  a result 
of  the  significant  figure  limitations  of  the  computer.  When  this 
happens  COMPE  is  too  small.  This  problem  can  normally  be  alleviated 
by  either  using  Precision  Focus  or  increasing  the  value  of  COMPE. 

2.  One  of  the  cable  segments  has  developed  zero  tension  or  has  gone 
slack.  In  this  case  the  array  is  statically  unstable.  This 
condition  also  results  from  no  change  in  imaginary  reaction.  The 
computer  program  for  the  tri-moored  array  has  been  equipped  to 


give  out  the  error  printout  for  the  benefit  of  the  user. 


In  Figure  18,  a logical  diagram, shoving  briefly  step-by-step  description  of 
hov  the  computer  program  proceeds  to  determine  the  equilibrium  configuration 
of  the  cable  system, is  presented. 
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Evaluate  the  weight 
forces  for  the  main 
array  and  assume  values 
for  FFX,  FFY,  FFZ 


Calculate  No-current 
configuration  under 
the  influence  of  forces 
FFX,  FFY,  and  FFZ  and 
gravity  forces 


Use  Subroutine  tie  leg  and, 
calculate  its  Equilibrium 
configuration  under  gravity 
forces.  Record  the  values 
of  RX1 , RY1  and  RZ1. 


Evaluate  hydrodynamic 
forces  under  assumed 
current  profile  and 
calculate  the  equili- 
brium configuration 
under  these  set  of  forces 
and  FFX,  FFY,  FFZ 


Use  Subroutine  tie  leg 
to  calculate  the  equili- 
brium configuration  of 
tieleg  under  the  action 
of  hydrodynamic  forces 


Test  for  force  balance 
criteria 


If  it  is  Satisfied 


If  it  is  not  Satisfied  - 
Use  Binary  Search  to 
calculate  new  values  of 
FFX,  FFY,  FFZ 


Store  the  No  - current 
(gravity)  equilibrium 
configuration 


If  it  is  not  satisfied, 
use  Binary  Search  to 
recalculate  new  values 
of  FFX,  FFY,  FFZ 


If  it  is  Satisfied 


Compare  the  No  -current 
configuration  with  the 
configuration  under  assumed 
current  profile 


Figure  lo:  Logic  Diagram  Showing  Step-by-Step  Description  of  Computer  Program 
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